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INTRODUCTION 


Tous les auteurs qui ont voulu étudier les*-noyaux atomiques 
contenant plus de quatre particules par une méthode plus précise 
que la méthode statistique ont employé la méthode du champ 
self-consistent. Toutefois ils n’ont utilisé cette méthode que d’une 
maniére approchée en construisant |’énergie du noyau au moyen 
des fonctions propres de Voscillateur harmonique. L’expression 
de l’énergie utilisée est souvent méme approximative en ce qui 
concerne les forces d’Heisenberg puisqu’on se contente de rem- 
placer les forces d’Heisenberg par des forces de Majorana d’am- 
plitude moitié (*), (), (*), @*), @), (*), @), @%), 7), 7), @), @), 
(7*), (), (*). 

Comme d’autre part la méthode du champ self-consistent n’a 
été développée que dans le cas des électrons avec action de Cou- 
lomb (*), (*), (4), (*), @°), 2°), 7), 78), 7%), (), (*), il était done né- 
cessaire de reprendre toute la théorie d’une fagon systématique 
en l’appliquant 4 un ensemble formé de deux sortes de particules 
distinctes entre lesquelles agissent des forces d’échange de diffé- 
rentes sortes. 

Ce travail se compose de deux parties. 

Dans la premiére partie qui est purement analytique, nous avons 
établiles systémes d’équations qui déterminent les fonctions d’onde 
individuelles et nous avons formé l’expression exacte de l’énergie. 

Nous avons exposé les principales propriétés des résultats ob- 
tenus, notamment la possibilité de mettre le premier membre de 
ces équations sous la forme condensée d’un opérateur unique pour 
toutes les fonctions d’onde d’une méme sorte de particules et 
utilisation de cette propriété dans un calcul de perturbations. 

Nous avons aussi montré comment se transforment les équations 
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et l’énergie lorsqu’on fait diverses hypothéses simplificatrices sur 
les fonctions d’onde individuelles. | 

Dans la deuxiéme partie qui est numérique nous avons expose 
une méthode complete d’intégration des équations du champ self- 
consistent. Dans cette méthode, qui est un perfectionnement de 
la méthode de Fock, nous avons conservé les Xx non diagonaux, 
que Fock supposait nuls, et nous nous en servons pour orthogo- 
naliser les fonctions obtenues a chaque étape du processus d’ap- 
proximations successives. Pour déterminer les fonctions propres 
et les valeurs propres des équations de Schrédinger nous utilisons 
la méthode de Milne (*) qui est trés rapide et trés précise et qui 
parait ignorée des techniciens du calcul numérique. 

Ces méthodes d’intégration sont appliquées au noyau d’Hélium. 

La valeur de l’énergie est améliorée de 8,5% par rapport a la 
valeur obtenue par Ja méthode approchée de Ritz-Heisenberg (*). 

La fonction d’onde finale, bien qu’ayant une forme analogue a 
celle de la fonction initiale, en différe, par endroit, de 7 %. 

La valeur de |’énergie obtenue en appliquant rigoureusement la 
méthode du champ self-consistent est de — 53,26 mc?. Pour rétablir 
Paccord avec la valeur réelle — 55 me? il suffit d’augmenter de 
1,16% le coefficient de la loi de force. 

Le peu d’importance relative de cette modification permet donc 
de supposer que la méthode du champ self-consistent est suscep- 
tible de donner des renseignements qualitativement bons surla 
structure des noyaux légers. 


ETUDE ANALYTIQUE DE LA METHODE DU CHAMP 
SELF-CONSISTENT APPLIQUEE AUX NOYAUX ATOMIQUES 


1, L’Hamiltonien des problémes nucléaires. Heisenberg (°), 
(**) a montré que les propriétés principales des noyaux atomiques 
sont explicables si on suppose ces noyaux formés de protons et 
de neutrons. 

La localisation de ces particules lourdes a l’intérieur des noyaux 
est concevable car les rayons nucléaires mesurés sont plus grands 
que la longueur d’onde de Compton relative 4 un proton au repos. 





h 
yt A ARIS 
A= a = 413 10-18 cm. 


Ces particules ne sont donc pas animées a l’intérieur du noyau 
de vitesses comparables a celle de la lumiére et il n’est pas utile 
d’employer la théorie de la relativité restreinte. 

Supposant le neutron formé par l’union profonde d’un proton 
et d’un électron, Heisenberg émit ’hypothése que l’action proton- 
neutron résulte d’une force d’échange analogue 4 celle qui assure 
la liaison homéopolaire dans la molécule H*. 

Dans la force d’Heisenberg les corpuscules en interaction 
échangent leur charge mais conservent leur spin. Le terme corres- 
pondant de l’Hamiltonien se composera done du produit d’une 
fonction de la distance des deux particules par un opérateur qui 
échangera 4 la fois leurs coordonnées d’espace et leurs coordonnées 
de spin. 

Majorana (?%) a proposé une interaction dans laquelle les corpus- 
cules échangent aussi bien leur spin que leur charge. 

L’opérateur correspondant n’agit done que sur les coordonnées 
d’espace. 

Toutefois pour justifier les énergies de liaison et les niveaux 
d’énergie des noyaux légers, ainsi que la dispersion des protons 
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et des neutrons, il a été trouvé nécessaire d’admettre que des 
forces analogues a celles que nous venons de décrire, existent éga- 
lement entre les particules de méme nature. 
Enfin Inglis (#2) a proposé, par raison de symétrie, de compléter 
Yopérateur d’interaction, en ajoutant aux forces déja connues : 
force classique : J(rav) +4 
force d’Heisenberg : J (rav)-P?%as 
force de Majorana: J(rap)- P7a0 


une action avec échange des coordonnées de spin seulement : 
J (Tas) + Pap 


Dans toutes ces expressions le symbole P représente l’opérateur 
de permutation ; indice supérieur indique la nature de la coor- 
donnée échangée: (x: coordonnée d’espace, «: coordonnée de 
spin) ; ies indices inférieurs caractérisent les particules sur les- 
quelles agit l’opérateur. 

Comme nous nous proposons d’appliquer systématiquement la 
méthode du champ self-consistent aux noyaux atomiques, nous 
adopterons ce type général d’interaction en particularisant la fone- 
tion radiale J(ras) selon le type d’interaction et la nature de la 
particule. 

En résumé, avec les types de force actuellement admis, l’ Hamil- 
tonien d’un noyau composé de n neutrons et p protons aura la 
forme : 





ena. Bee ee ea Ai 
7 8x2m, sd ‘ 8x2m, sand : 
t=4 i=n+1 


hs 3 [Gta(riz) + GPn(riz)Pix + GPn(riz)P%e + G4, (riz) P?%ix] 
i>k=1 
n+p 
- 
+ WS (G(r) + G2,(re)P ie + G*,(rit)P%R + G4a(riz)P* ie] 
t>k=n+1 
ia ae 
+ > »} [J*(rit) + J?(riz)Ptiz + I8(rix)P Si, J*(rix)P2>44). 
t=1k=n+1 


Dans cette formule, m, est la masse du neutron et m, la masse 
du proton. 


Les coordonnées des neutrons sont numérotées de 1 a n et les 
coordonnées des protons den +1an + n: 
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Nous désignons par la notation abrégée 


n n 


n 
la double sommation » 
i>k=1 k=1i=k+1 


Entre les différentes sortes de sommation existent les relations 
suivantes : 


n n n nm n 1 n nm 
FR ly ela oN 

. ~ y, aut 
t>k=1 k=1i=k+1 = i=1 k=i+1 i=1 k=1 
ke% 


| 
i=1k=1 i=k=1 


2. Principe des méthodes d’approximation (°*), (*), (14), (1°). —Soit 
un opérateur hermitique H. 

Nous nous proposons de déterminer les fonctions VY normalisées 
qui rendent minimum l’intégrale 


fwsttve: 


Pour cela nous devons résoudre le probléme de variation 


I 
Le) 
C | 
M 
Ve 
t4= 


(2) 5 [ wertwas = 
avec la condition 


(3) | wud: = 1. 


Selon la méthode habituelle du calcul des variations nous devons 


annuller 
5 { wentwas — 2 [ vere: 
ce qui donne |’équation : 
(4) fare — E)Wd: + fora — E)é¥dt = 0 


Remarquons que E étant une constante réelle, H — E est un 
opérateur hermitique. 
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Mais d’aprés une propriété connue des opérateurs hermitiques 
fovea — E)Wd =f oven — E)*0*dz 
et Péquation (4) devient : 
ii aw*(H — E)Wde + if éw(H — E)*w*dz = 0. 


On voit que le deuxiéme terme du premier membre est |’ima- 
ginaire conjugué du premier terme, et |’équation (4) peut s’écrire : 


(5) partie réelle de [ f. éw*(H — Eyvat | 0. 


Cette équation (5) sera a fortiori satisfaite si on pose simple- 
ment : 


(6) fovea — E)Wd: = 0. 
Si on n’impose aucune condition restrictive 4 la fonction ¥, 
la solution de (6) s’obtient en résolvant l’équation 


(7) (H — E)v =0 


qui n’est autre que l’équation de Schrédinger. 

La résolution de l’équation (7) donne la solution rigoureuse du 
probléme physique correspondant A Popérateur H. 

Mais il arrive fréquemment que l’équation (7) soit trop compli- 
quée, et on se contente d’une solution approximative. 

Pour montrer sur quelques exemples comment on cherche une 
solution approximative, supposons que l’opérateur H porte sur 
deux variables z,, et z,, et cherchons une fonction dela forme : 


V(x, %2) = (2): }(X2). 
qui satisfasse a l’équation (6). 
Cette forme de fonction d’onde se justifie notamment si Popé- - 


rateur H est symétrique en XZ, et Lz. 
L’équation (6) devient : 


(8) { ey (ay) -Y* (a) (H — E)¥(cr) (x) dard 


+ [ere Wet Pied E)¥(2,) -V(,)dayde, = 0 
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A cause de la symétrie de H, la deuxiéme intégrale est identique 
a la premiére. L’équation (8) se réduit done a: 


ah by*(a) «Y*(a’)-(H — E)¥(2)-4(a')dada! = 0 


et l’équation d’Euler de ce probléme de variation s’obtient en 
annulant le coefficient de d)*(z) : 


9) iE Y*(e’)-(H — E)¥2')-4(a)da’ = 0. 


L’intégration éliminant la variable zx’, l’équation (9) ne porte 
plus que sur une seule variable, et la difficulté est évidemment 
diminuée. 

On peut aussi chercher une fonction de la forme : 


W (ay, Xq) = Y(X, — 2p). 


Ce serait le cas si Vopérateur H comportait des termes relatifs 
a la distance des particules M, et M, de coordonnées 2, et 2;. 
Dans ce cas il serait indiqué de faire le changement de variables 


Ly —t =U 

yi+%my=-? 

Inde , ed D(z, Lo) ais 

Pon: dx,-dx, = (eR ; ) du do = 
L’équation (6) devient : 


4 
5 du dye. 


f eeu) cH — E)-(u)- 5 du- do = (0) 


on en déduit 
[ (1 = B)-4w)-ao =0 


et. dans ce cas encore l’équation du probleme ne porte plus que 
sur la variable u. 

Considérons le cas ou (H — E) opére sur les coordonnées po- 
laires (p, 9, y) d’une particule. Pour simplifier, et si la nature du 
probléme physique s’y préte, on peut décomposer la fonction 
d’onde en un produit d’une fonction du rayon par une fonction 


des angles polaires 
¥(9, 9, 9) = fle) ¥(9, #) 
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la fonction f() étant l’inconnue a déterminer, et la fonction 
Y(9, ~) étant une fonction sphérique connue. 
La relation (6) devient : 


J torr, ¢)-(H — E)-f(2)- (0, ¢)-p%dp sin 6-dd-dp = 0 
et ’équation du probléme est : 
J Y*(6,)-(H — E)f(e)-¥(0, ¢) sin 6-d0-dy = 0 


c’est encore une équation & une seule variable p. 

Il est évident que plus on précise la fonction Y, c’est-a-dire 
plus on diminue la part d’inconnu qu’elle contient, plus on facilite 
la résolution du probléme, mais plus on s’éloigne de la solution 
rigoureuse. 

A la limite extréme de ce procédé on trouve la méthode de 
variation de Ritz d’aprés laquelle la fonction ¥ est complétement 
fixée dans sa forme mais contient simplement des paramétres dont 
on détermine la valeur pour satisfaire aux conditions (2) et (3). 
7 D’aprés (7) on voit que 


if v*HWd: = E = énergie du systéme si H en est l’Hamiltonien. 


La condition (2) s’écrit done 
cE = 0. 


Si la fonction Y qu’on s’est donnée a priori dépend de p+ 1 
parametres : a, o%,... %p,,, Pun de ces paramétres, par exemple 
ap,, est déterminée par (3) et les p paramétres restants servent 
& minimiser l’énergie. Ils sont done déterminés par les p relations 


aE ey gees 


(10) =, = 9, 
0% 0% op 

Comme exemple d’application de cette méthode considérons le 
cas ol on connait une suite compléte de fonctions orthogonales 
portant sur les mémes variables que la fonction inconnue Wet 
satisfaisant aux conditions aux limites et aux conditions de conti- 
nuité imposées 4 ¥. On peut se proposer de développer en série 
Y suivant ces fonctions. Si on s’arréte au (p + 4)*™* terme on 
a précisément la fonction approchée ¥ de la méthode de Ritz, car 
les (p + 1) coefficients inconnus du développement sont les p + 1 
parametres 4 déterminer par les équations (10). 
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3. Méthode du champ self-consistent (*), (°), (#), (45), (%). — La 
méthode du champ self-consistent est une méthode d’approxima- 
tion destinée a simplifier le traitement des systémes de corpus- 
cules en interaction. 

On sait que si plusieurs systeémes sont indépendants l’ Hamil- 
tonien de l’ensemble de ces systémes est formé par la somme des 
Hamiltoniens de chacun d’eux et la fonction d’onde globale est 
le produit des fonctions d’onde partielles. 

On voit en effet qu'une fonction de la forme 


(11) W(Xq, Vqy°* +n) == Vy(2,) XK Vy(Xq)+ ++ X Vn(In) 
ne satisfait a l’équation : 


(12) HW = (Hy(2,) + H,(2_) ++++-+ Halaa))¥ = EW 
que si ¥,, ¥,,.. Yn et E sont déterminés par : 

H,¥, = E,", 

HV, = E,V, 


Ha¥n = En¥n 
E = E, + E, +---+ En. 

Ce résultat ne se conserve pas si les systémes considérés sont 
des corcpusules identiques. Dans ce cas la fonction d’onde globale 
doit étre antisymétrique pour satisfaire au principe de Pauli, et 
il faut la constituer par une combinaison linéaire antisymétrique 
de termes analogues aux précédents. On sait qu’une telle combi- 
naison est identique au développement d’un déterminant 


(43) v= ; 
Wy (an) + + + Vn(an) 








Chacun des termes de ce déterminant satisfait bien a l’équation 
(12) car les Hamiltoniens des différents corpuscules sont identiques 
et les fonctions Wi sont fonctions propres d’un méme opérateur. 

Pour cette méme raison on peut aussi supposer les fonctions T's 
orthogonales et normalisées. Dans ce cas on normalise la fonction 


tee ; : 1 
globale Y en multipliant Je déterminant par oe: 
nm 


La supposition essentielle qui caractérise la méthode du champ 
self-consistent est que les corpuscules sont indépendants. 

Par conséquent la fonction d’onde de ensemble doit avoir la 
forme (13). 
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Parmi les fonctions de cette forme, celle qui s’éloigne le moins 
de la solution exacte est déterminée par Ja méthode de varia- 
tion. 

“On aboutit ainsi A un systeme d’équations intégro-différentielles 
qui détermine les fonctions ¥%. 

Ce systéme d’équations a une forme telle qu’on peut linter- 
préter en disant que chaque corpuscule se meut indépendamment 
des autres dans le champ moyen de Il’ensemble. 


4. Etablissement des équations des fonctions d’onde partielles 
(*), (®). — Considérons un noyau constitué par m neutrons et p 


protons. 
La fonction d’onde du systéme a la forme 


(14) y= VY x WV, 
avec 
1 (74, %) +++ Yn(2y, 2,) 
(15) De eet fF OR ee 
Vn pfs Sisktth © sie a erage oie eels 
Y (a, Tn) ++ Un(Xn, Sn) 
1 Ya(Tn1, Ons1) ee Up(In41, n41) 
(16) We tena | 0b See oa : 


Vr(Lnsp, ony p) ae ad $r(Lnip, Cnsp) 


Nous désignons par a l’ensemble des coordonnées d’espace x, 
yi, z eb par a la coordonnée de spin de la jime particule. 

Les coordonnées des neutrons sont numérotées de 1 a 7 et celles 
des protons de: n+ 41an + p. 

La fonction ¥, et les fonctions Yi caractérisées par un indice 
minuscule sont relatives aux neutrons. 

La fonction ¥, et les fonctions Y, & indice majuscule se rap- 
portent aux protons. 

Les fonctions ,, Ye,-. et dn peuvent étre Supposées toutes 
orthogonales et normalisées, car elles peuvent étre remplacées 
par n fonctions orthogonales obtenues par substitution linéaire 
puisquwune substitution linéaire quelconque parmi ces n fonc- 
tions a pour effet de multiplier le déterminant par un facteur 
constant. 


On peut supposer aussi les fonctions orthogonales et norma- 
lisées. 
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L’orthogonalité et la normalité des fonctions d’onde partielles 
s’exprime par la relation 


#( a). Ache Poke 1sit=k 
> fe (x, )+ve(a, ) da f= Osiixk. 


= 3 


Appliquons a la fonction ¥ la méthode de variation définie 
psécédemment, en adoptant pour Hamiltonien du systéme l’ex- 
pression suivante : 


pat Mace 
Be kag i! Bate, 
i=1 i=n-+1 


n 


+ >) [Ga, x) + Gai, x)P% + Gai, rx)P% 


t>k=1 
+ G4(axi, xe)P ix ] 
n n+p 
+> Dd [Pes 2) + Me, PH + Me, m)P3, 
i=1k=n+1 


+ J4(xi, ae) Pi | 


Pour ne pas alourdir l’écriture nous ne représentons pas les opé- 
rateurs d’interaction entre protons car ils sont tout a fait ana- 
logues aux opérateurs d’interaction entre neutrons. 

Quand nous voudrons rétablir la symétrie des résultats nous 
désignerons par Gn', Gr?, Gn, Gat les opérateurs d’interaction entre 
neutrons et par Gp!, Gp?, G,*, Gp*les opérateurs analogues entre pro- 
tons. 

Dans le méme but de simplification posons : 


h? 
Sat el: : A; = Hi 
h2 
es ond 2 
— or ‘ Ai — H; . 


On peut décomposer I’Hamiltonien de la fagon suivante : 


H=L,+L,4+L1,+ ly 
avec: 
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n+p 
| 
y He 
i=n-+1 
n 
L; = )) (Ga, a) +--+] 
t>k=1 
n n+p 


> [J*(xi, Xe) +--+]. 


i=1k=n-+1 


Ly 


Les équations qui définissent les fonctions d’onde partielles sont 
les équations d’Euler du probléme de variation 


Ds fever — E)¥da: = 0. 
O1,°**Onip 
Orie 
owe = OW *.W,* 4 W*.dW,*, 
La variation totale se décompose donc en deux variations par- 
tielles 0,1 et 0,1 


6,1 = >; f aewn(H — E)¥,-W,-dt,-dz, 


Ty," * "Snip 
Oy,°° ‘onip 


Caleul de o,f. 


(17) 61 = » f eu Hyh Ly — B+ Ly Ly + Ly) hy Wye 


Sy,°° *Cnip 


>» fmt, — E) + B, + A,L; + D,]-Wydr,. 


I 


Nous avons posé, pour abréger 


es p> fees, A,= > fv tvy dt: 


c ee 
1s°°°On Sn+19** * Snip 


Ba 2 aS ‘Ly-¥,-dy,B,= J *LW adr: 
Sn+15°* * nip 
os y [ ¥y*- Lye, 
Ar: ad 


= Wr. *, S on, 
: > if Fy*-L,V,-dy, Dy = > aE V2" LV dt, 
1°" "Sn Sn+19"** Sn+p 


dt, => dz,,: a ‘AXn, dt, => ditn,1: Mg *Atnsp- 
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Calcul des termes A. 


Ay= »~ f wt ry dey 
F,°* "Snip 
On doit prendre dans ¥’, et ¥',* des termes provenant des mémes 
permutations des indices, car deux termes provenant de permu- 


tations différant par exemple par l’échange des indices r et s 
donneront naissance a l’intégrale 


{ 
ratte Lf Yr (ai, 64) -Yo¥(ae, on) Yalri, 2%) Yelm, cx)dae, dre 


95,57= + 


l__ 


~ 


4 ‘ — > 
= ib f He les ander x S Jf Ytlen, 00) tela, o4)- dae = 0, 
O5 Pe 


a cause de l’orthogonalité des fonctions ¥. 

D’autre part lintégrale portant sur le produit d’un terme de 
'’,* par le terme correspondant de ¥, est égale a l’unité car les 
fonctions ‘’ sont normalisées. Comme il y a n! termes dans ,, 
on voit que: 





19 A, =nix Pea e os 
ia : Vn! Vn! 
de méme 

(20) A ae 


Calcul des termes B. 


avec 


L, = » He. 
i=1 

Comme l’opérateur L, n’agit que sur les coordonnées d’une seule 

particule a la fois il faut prendre dans ’,* et Y’, des termes corres- 
pondant 4 la méme permutation. 


Un terme quelconque de ¥',* donne naissance aux intégrales 


n 
a » ff ps e¥ (tn, Fm)+* = Ye (ey o)* Het bem, Im) ++ Ye(Ze, e)Ea, « « - 


-+A2m:+-OxLE-* -dIn. 


M. MATRICON 2 
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Une de ces intégrales s’écrit par exemple : 


> fw (Xm; om)*H m 1) (Lm; 5m) dL = 


mits 
mae 


2D fete ote, oe, 


A 
Pies 


car le nom de la variable d’intégration n’a pas d’importance. 

Pour terminer le calcul nous ferons la remarque essentielle sui- 
vante qui servira constamment par la suite : 

Prendre successivement dans un produit de n fonctions de n 
variables de la forme 

Yr(Tm)- + -Ys(Tn) + - > Ye(ai), 
toutes les valeurs des indices des variables revient 4 prendre 
successivement toutes les valeurs des indices des fonctions. 

Par conséquent dans l’opérateur L, la sommation effectuée sur 
Pindice de la variable peut étre remplacée par une sommation 
effectuée sur lindice de la fonction. 

Le terme considéré de ¥',* donne donc naissance a |’expression 


n 


1 
ni’ 2a Xf ee 9 Wve, oa 


= +5 


On voit que Pindividualité du terme considéré de ¥,* n’inter- 
vient pas dans le résultat ; par conséquent les n! termes de SP 
donneront une contribution identique et au total on obtient 


n 


B= DL fete. 9-1-we, oe 


om 5 


Par analogie 


Pp 
B= >\- 1 J (x, c)- H?. *Uy(2, o)dx. 


Une intégrale de B, n’est autre chose qu’un élément de la ma- 
trice représentant l’opérateur H! dans le systeme des fonctions ¢. 
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Introduisons donc les notations 


Hi, = x fe (x, )+H?-dy(a, o)dx 


o=45 


Hix = : Jes )H*y,,( (x, s)dx, 


e=+5 


ce qui permet d’écrire 


ie = > Hi 
i=1 


(24) : 
. 
| B, = > Hj. 
I=1 
Calcul de C. 
Gis “ if W,*L,V dz, 
TytsSy 
avec 


Dd (Gee, re) + Gai, 24)P%, + Gai, m)P%, + G4(ai, m)P72]. 
i>k=1 
Pour effectuer le calcul décomposons l’opérateur L, en opéra- 
teurs partiels 


n 


[Ln = >) Ges, ms) 
i>k=1 


Liss = 2 G?(2;, %)P 
ike 


n 
L33 = D3 G(2i, te) Pix 
i>k=1 


Lig, os G(x, x)P 
i>k=1 


Le terme C se décompose d’une fagon analogue 
C=C, + C, + C;, + Cy 


J n 
Gis yy fu SY Gia, 22) ¥ de. 
34° +On t>k=1 


(22) ( 
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Les termes provenant d’une méme permutation des indices dans 
y, et ¥,* donnent une contribution au terme C, qui a pour valeur 


n . 
VOY | te) ete, 2), 2!) h(x, 9) Yale’, o'derde’. 
i>k=1¢,0' 
Si on considére maintenant un terme de ¥,* et celui de '¥, qui 
s’en déduit par l’échange des particules i et k, on obtiendra une 
contribution au terme C, de la forme : 


Jel) = > f ee 04) Vs" (Xx, oz) -G1(24, rx) 
S5,5k 


Ys(2i, %) -U(Xk, cx) day: dry. 


Le signe moins provient de ce que les termes de 'V,* et de ¥, 
correspondent a des permutations de parité différente. En faisant 
la somme des expressions analogues pour tous les termes de ¥’,* 
et pour toutes les valeurs de i et k, on trouve 


Sha 
i>k=14¢,0' 
Ye(2, s) ; hi(z', a’) -dxdzx’. 


Posons : 
Sf v.*. 0) aXe’, 0") -GUa, 2')-t(a, 0) dale’, o')dxde’ = 
(tz G! | t 
ay than ae: 


La dépendance entre la position des indices des fonctions et la 
nature des variables est expliquée par le schéma ci-dessous : 





y* Y 
SS —— a, ei eee oe 
' t i i } H 
I 1 1 I ! i 
RN Tian fee a ! 1 rat 
| } Phare ie Poa I 
Rina, Breen cts, 
TEMG ELLE hy as Sy 1 ' { 
Si 
o tot fous o 
H i 
SSN SE HS ESE penne | GE | met 
a’ k { ; 
ae : 
i { H 
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On a done 
Se leet alo | ll 
On trouve facilement de la méme facon 


om SS (ale fal ll 
an BS ftel 1A lal* Lat] 


oe a a a 








Calcul des termes D. 


D, = Sy f vtbends 


Oy+++on 


Décomposons L, comme nous avons décomposé L,. 


| n n+p 
%. 2, 
La =» >» J*(xi, Xx) 
i=1k=n+41 
n n+p 
a4 x : 
Lae >, dy It; m)P ip 
i=1k=n-+1 
(23) n n+p 
I45 =>) YS Iai, we) Pix 
i=1k=n+1 
no n+p 
Digs =» p> J4*(ai, rx)P i; - 
i=1k=n+1 


Remarquions que l’opérateur L, n’opére a la fois que sur une 
vet 2 
seule variable de W,, nous ne devrons donc considérer dans ¥’, 
et ¥, que des termes correspondant aux mémes permutations. 
La décomposition de l’opérateur L, divise D, en quatre parties. 


D, = Di; Tr Die + Di3 no Dia 
no n+p 


=, me 
Pa= DY ford Y Ben avid, 
Tyson i=ik=n+1 
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Un terme quelconque de ’,* donne naissance a : 
n-+-p n 
f ¥ * 1 
— YY] dS vee, Ie, aa)yila, )dz. 
1 


n! 
k=n+1 o ix 


Pour tous les termes de ¥,* on obtient donc 


n+p 


Du = » aS Yi * (x, s)-J4(x, Le) Ya(Z, s)dz. 
k=n+1 o t=4 


On ne peut pas calculer immédiatement les autres parties de 
D, 4 cause des opérateurs de permutation de coordonnées. 
La seule partie du terme D, calculable actuellement est : 


n at 
Da = J & Ya¥(x, 2)-I'(x, ai) velar, o)de. 


i=1o 


Ces calculs préliminaires étant terminés nous pouvons récrire 
la valeur de 0,1 en tenant compte de (17), (18), (19), (20), (22). 


6,1 = (Bg — E) » faye 
Tyo on 


+» af oW,*L, dz, 


Oy ° on, 
a yy sj oW*L5¥ dt + » i 61 *D.V dz. 
Oy" ++ On Cys oy 


Calcul du premier terme de 3,1. 


> fet ds, 


1 . "on 


3,1, = (B, — E) 
3 
Considérons un terme de ¥’,*. Pour faire varier ce terme qui est 
un produit de fonctions d’onde partielles, i] faut prendre successi- 


vement la variation de chaque terme du produit, c’est-a-dire faire 
sur lindice i une sommation de la forme 


n 
Mey (a). 
i=1 
L’opérateur 0 est donc équivalent a un opérateur agissant sur 
une seule variable a la fois. I] faut donc prendre dans ¥',* et ¥ 
1 
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des termes correspondant a la méme permutation. En tenant 
compte du fait que les n! termes de V’,* donnent le méme résultat 
on obtient : 

n 
‘~) “ 
ah = Br —B)D f Yara, 9)-tle ade 
¢ = | 


Calcul du deuxiéme terme de 6,1. 


Soon my forth 


Oy on 
‘ n 
7 “| 
te »> four M3 HW dx. 
Oy++*Gn er 


Considérons d’abord dans ¥',* et Y, deux termes provenant 
d’une méme permutation. 

Les opérateurs 3 et H! peuvent agir soit sur la méme variable, 
soit sur des variables diférentes, ce qui donne des termes de la 
forme 

Bta* (ar, Sr) H+ (ar) Yi\ tr, Or) der 
et 
SYi* (Lr, Sr) -Ug*(Xm Fm): Yi(r, Fr) - Hy (2m)- Ye(Lmy Fm)ALrALm. 

En faisant les sommations sur i et sur m, ou, ce qui est équiva- 

lent, sur i et sur k, et en tenant compte des n! termes de V',* 


on trouve 


{ n 
\ Df Deere) eH de 
¢ i=1 


| +> i Seutiese) SL Yate, eV Wale edede. 
¢,0'* i=1 


mf<i=1 


(24) 


Considérons dans ¥,* et yy’, deux termes qui différent l’un de 
autre par l’échange de deux variables. Pour que l’intégrale por- 
ermes ne soit pas nulle il faut que l’opérateur 


tant sur ces deux t 
opérateur H, sur l’autre. On a donc 


) agisse sur une variable et I’ 
un terme de la forme 


_ a byi* (Xr; Gr) t Ux(Xr, Gr) f v,* (2s, Gs) H1(2s) : Vi(Ls, Ss)AX,y -d2;. 


: Gr,Ss 
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En faisant les sommations habituelles on obtient la troisiéme 
contribution a d1,, 


(25) —) if Sevi*(a, 2) M3 Ye(X, o)-be*(2', o')-H1(a')-di(2’, o’)dxda’. 
o, a! i=1 kAi=1 
Dans (24) et (25) remplacons 
> Sy, fix par ¥| Se — fo | 
t=1 k~i=1 i=1Lfk=1 


n 


on voit que les termes en y, fi provenant de (24) compensent 
i=1 


exactement ceux provenant de (25) et il reste finalement : 


‘1, =>) | Ser *) | | w + Yi [ote e) —) Htada(a, 2) | dr 


k=1 c—4 


Calcul du troisiéme terme de 0,1. 
61, =™ f ov, *L,V dz). 
Bisson 


Ce terme se subdivise en quatre a cause de la décomposition 
(22) de L,. 


n 
BT ote b 4 
Phi, ow,* : G(2xi, x4) Vd. 
Ty s*+ on iSke1 


Ecrivons l’expression a calculer sous la forme équivalente 


n 
1 & 
5 * > iy ‘ 
S11 A Nk eh 3 Gl(ai, ag) Vdzy. 
Oy°* +n i~Ak=1 
Considérons un terme déterminé de ¥,*. Nous savons que tous 
les n! termes de ’,* donneront le méme résultat. Nous multi- 


plierons donc tout de suite nos calculs par n! ce qui a pour effet 
: : i vhs 1 
de faire disparaitre le coefficient as 
n?i 
a) Prenons le terme de VY, qui correspond a la méme permuta- 
tion que le terme de \* considéré. 


1° L’opératieur 3 et Popérateur G agissent sur une méme par- 
ticule. 
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Les termes correspondant de 01,3, sont 


| 


| =. 
li x by,* (x, s) r Ys*(2', 3’) Gl(z, x’) : Ys(2", o’) ; (x, s)dxde' 


| 


xfs D a*(a, 9) dole, ¢) G(x, x')- Mure, 5’) U(x’, o’)daxdz’. 
o,¢ ence ar 
aap eS = OYs* (x, ¢) + Ye(x, o)-G1(ax, x')hs*(2', o')-Y,(x' o')dadz’. 


2° L’opérateur 3 et Vopérateur G agissent sur des particules 


>! > Sy ie tira hi 2) Yala, 0) deed’ 
areal 


(26) 


différentes. 
Un terme s’écrit 
2 _f entte )-e*(a', 0!) -bu*(x", o"")GH a’, 2”) -e(x, 2) + ¥,(2", 0’) 
oe B bu(x”, o”\dada'dx". 
=“ 6,* (zr, c) atl Gy he t )%s(x, s)dx 
c 


Pour obtenir tous les termes de ce genre il faut effectuer sur 


Vintégrale précédente les sommations 
n 


my Xs: 


n 
N 
eed ya 
4s=1 utiss=1 


No — 
a sr 


On trouve aisément que ces sommations sont équivalentes aux 


suivantes : 
n n n n n n n n 
4 N te NV be . 1 . 5 Vv 
27). = / —F Sea <a / + ag" 
! ia 2 s=1 t=u=1 2 t=4 i opp 2 u=1 Arve s=t=u=1 


s,t,u=1 
L’ expression sommeée étant symétrique en ¢ et uw certaines somma- 


? 
tions telles que les suivantes sont identiques 
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et finalement la contribution a 41,3, des termes considérés a pour 


expression 
n 
5m if Sane X, 5) 








pa | te | GI | Ey vs(x, o)dx 
s=1 t,u=1 
3) Yar 2) >») {+ | Gi t Wha E 
t=1 
— Df Sar (x, s) >> Ve | Gl i. Us(x, o)dax 
s=1 t=1 


+ Ef dares (x, s) 


+7 | Gi tr me b.(x, o)dx. 


b) Prenons dans ’, un terme qui différe du terme considéré 
dans ¥,* par ’échange de deux particules 


Il faudra donc mettre le signe moins devant tous les résultats, 
ce que nous ne ferons qu’a la fin 


1° L’opérateur 6 et lopérateur G agissent sur la méme particule 


On déduit immédiatement de (25) par échange des coordonnées 
dans les fonctions ¥. 


Sf yews Z, 5) +(x’, o')-G1(a, x’) Son x',5')-U,(ax, ¢)dadz’ 
s=1 


s=1 


-J Sapte, 9) Ys(z', o')- G(x, x')-b*(2', 0’) -Ye(x, o)dadx 


2° L’opérateur 0 et Popérateur G agissent sur des particules 
différentes. 


Il existe trois sortes de termes 


2 


pas O}e*(x, 0) -Ye¥(a!, 0!) + hu*(x", 0") Gar’, a") hex, 5) -Y,(2", o”)- 
er ne Oe*(2, 5) Ye (z', 0!) -by*(a", o”)-GUa', 2") -de(a', 0’) 


eo f euttao 


Yu(x', o')dada'dz’. 





(x,0)dx 





Ye(a, ¢)- Yu(x", o")dadz'dz". 
t | Gl | BN | Yle, s)dx 
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> [ante 3) 7 be*(2’, a’) . hu*(2", o")Gi(2', x") Ye(x", 3") i Y.(2, 3) r 
aes by(a’. o’)dada' da’ 


=ps i ay,*(a, 2) | t | G | uw | U(x, o)dx. 


Il suffit maintenant d’effectuer les sommations (27) sur ces 
termes. Nous n’écrirons pas les résultats qui sont trop encombrants. 

c) Prenons dans WY, un (terme qui différe du terme considéré 
dans ¥',* par l’échange de trois particules. 

On affectera le signe plus aux résultats car les permutations 
des deux termes sont de méme parité. 

On obtient les deux sortes de termes : 


fetta 3) J U*(2', 3’) y hu*(2", o")Gi(z', x") -Y(2', 3’) z Y(x", 3") 
v(x, )dxdx'd' 


‘ 
,c/,o" 


Efren 


fave (z, s) . Ye* (x’, a’) . bu*(2", 3") , Gl(z', x’) i Y,(2', 0”) : Ye(2, 8) 
nal vu(a’, e')dar- da! - da! 


-> if outta) {t |G) a | ula, dee 


Finalement, en désignant par ») ensemble des sommations 


-by(x, 0)dx. 





ay iu 
(27) on obtient : 
(28) dl, = 
ih Sabet (oo) tle #) GH, 2) als 9) He ee 
oc t,s=1 


Sal: Sayer(a, a)-Ys*(2', s') Giz, x!) Ys(, 9) ° Y,(x', 0 dada’ 
3,o! s=1 

LS Sere 9) he 2) GG 2) 0) He a)diede’ 
oa’ t,s=1 


+ yl bYs* (2, ¢) ; hs* (z'; a’) ; Gi(x, x’) r (2, 3’) , Ys(2, s)daxdz' 
3! s=1 


5, 
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Poa Mey *( (2) [{* ioe 
u 


ene 





t s(x, 6) 
u 
as ie te GI | hea! 


et ha ina 











not | t |G|t [ela s) 
u $ 
tue rhea 
oh | é G! a v(x, 2 ba 








En appliquant la méme méthode on obtient les autres contribu- 
tions 4 0,1, 


(29) BIgo 


sf So, *(x, ¢)-Ye*(2x', o')- G(x, a’) -U4(z', 0) - U(x, o’)daxdz' 


—, 


pal Yawe, c)+¥eR(x!, 0") -G2(x, 2) Ye(a,¢) dal, «dard! 


vi DED? y bhs* (x, o)-1,*(x', o')-G 2(x, x’) «s(x, 0’) U(x’, c)dxrda! 


let 


+e x ste *(a, 2)-Ys*(a', o') -G2(x, 2’) U(x, 9°) -b5(x’, )dadx’ 


+L J Saree lee 


S,t, u, 


(G2 





ie Us(z, c) 





af 


mean 
u 


— jt jo 





G? | u }s(z, 2) 
t 


G? 





f U(r, o) 








l bu(x, 3) 
s 


+ | ‘ G? | . bu(2, c) 





“+ . | G? | mh U(x, ?) dx. 
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®1133 = 


sh Nove (a, ¢) = be*(2', a’) G(x, 2) : d(x, 3’) -U,(2', «)dadz' 


t,. s=1 


Sai Sey (w, 9) YX’, o')- GPx, a!) Yale, 2!) Yala", 2)dndar 

c o! s=1 

-yf Miye*(e, a): U;*(a', o’) -G3(z, rye Us(a", 2) -Ye(a, a \dxdsr' 
a,c s,t=1 


n 
+Ef Sey.*(x, a): s*(2', a’) : G(x, a’) ; U(x 3) Ys(z, a')drdx’ 
oo! t=1 




















= { Laven (|| | bose 
& _- t G ie | Ys(2, ¢) 
o G oH 3) 
u s§ 
rae eaten s) 
+E 1e ig AAS 
Ss 
aE G3 a | le fr 
u 








Oly, = \ 
Bi Dive" (x, ¢) . Y* (z', ’) : GA(x, x')¥s(x', o’) ; U(x, o) dada’ 
o,0’'% t,s=1 
n 
->{ Saye*(e, 9) etl 2) Ga, a! ola’, 2!) tale, e)aenda! 
¢,o s=1 
n 
yf Syste, 2) dt al 0!) GMa a!) Yolo, 2) Ha! o deed’ 
co,c'% t,s=1 


Laps if Say (x, 2) Yala’, 2!) GA, 2!) Hele, 0) -Ye(2', s’)dardar’ 
3,3! s=1 


+> up DPv.t(25 2) | ” | G* 


bate Gt 





u [He 
t 
t U,(2, 2) 
u 
Abate 
— i G4 |s po) 
u t 
+(t Beh \ Yu(x, 9) 
u Ss) 
Suk G4 | s i 
u 











u 
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pees 

Dans les expressions (28), (29), (30), (31) que nous venons d’écrire, 

le deuxiéme et le quatriéme terme s’annulent. D’autre part, en 

explicitant les sommations d’aprés (27), et en réduisant les termes 

semblables on obtient les variations sous une forme plus condensée. 
a1, = 


n 


nai Soye*(e, g)- po pe*(x', 3’) 7 pa’, 3’) 4 G(x, z') : Ys(2, 3) 
3, a! s=1 t=1 U(t, eZ Gia, z')- bs(2, a’) 
-G? 


~_— 


x, z')+Y,(x', c) : 
x, 2’) -bs(x, 0’) dxdx 


z eo: Yel, a’) 


nn 
* %* 
om FM; 
k= tat 
a a, 9 
— — ~~" 
ce 
a Ks 
om <A Poplar caer 
Tae 
wo wo 
o~ om os 
oe 

S 

= 

< 

a 

=— 

& 

2 

etl 


+ Y*(x', 0) -Yu(x, 5) Ge, x')- RC ‘ 
pao Uy*(2', a’) d(x’, a’). G4 (x, z')- -b4(Z, 3) 


Sf Spend 
s s=1 


n 















































ABO Me as aoe Henn) 
sao dd 

i belahelbak S: mar eps ia G 2) me) 
—{t, Gi igus G? tle GS tre Ae G*|u i" 


























Calcul du quatriéme terme de 5,1. 


a1, = 3 { bY FW ,*L WW dz, «dt, 
Tye Onin 
La décomposition (23) de L, en quatre termes a pour effet de 


diviser cette variation en quatre parties qui ont les valeurs sui- 
vantes : 


n n+p 
. . 9 = : 
Oly = ~ fomewst), ° S: J1(xi, tz) V+ Wd, - Te 
Ty" *ontp i=1 k=n+1 


L’opérateur qui agit sur Y, étant Popérateur J*(ai, 2x) qui ne 
modifie qu’une variable zz a la fois, il faut prendre dans ¥’,* et 
, des termes correspondant A une méme permutation. 

a) Considérons d’abord dans ¥,* et ¥, des termes correspondant 
& une méme permutation. 
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1° Les opérateurs 3 et J, agissent sur la méme variable dans 
WY * et Y;. 


La contribution correspondante s’écrit : 


n p 
oA) SD syste, 2) tact a’, Ma, 2°) -Yele, 9) ule’, 9!) dda’ 


s=1M=1 


2° Les opérateurs 3 et J‘ agissent sur des variables différentes 
dans ¥,* et ,, ce qui donne : 


n 


Sy Sra (x, 2) a(x, 2)- S Yek(a’, o') we’, 0") - Ma, 2” 


a, o! vase ttAs=1 


P 
Sate’, of) ula’, o) ded! da! 
M=1 


n P 
= 2 si ao >? e*(ar, 0) =Ys(2, 2) -Ye*(a', 0’) Yuk (2", 3”) -I*(2', 2") 
g, a” vg? 71 M=1 
Yel’, o')-ru(2", *)dada'dar” 
P 
Ri » ps D> du-*(e, a) +Ye(x, 2) > e*(x', 3’) -e*(x", 0”) - IM’, 2") 
s, 3’, s=1M=1 
(2, a’). Vy(x" 3" )dxdx'dx" 
b) Considérons maintenant dans ¥',* et des termes correspon- 


dant a des permutations différant par l’échange de deux particules. 
On obtient immédiatement 


P 
Ve >» Ai >, Sev*(2, c)-h,(2', a’) UF (2', a’) du*(2", a") 


s,t=1M=1 
Ji(x', x) . H(z, ) Yu(", a" )dadx'dx" 
n Pp 
> ips dav." (2, s) “Y,(2', ’) “hs*(2', 3’) “u*(2", s") 
3,0/,07" s=1M=1 
F(z’, x’) Xs bs(2, s) ‘ Yaz", o")dadx'dx". 
La deuxiéme partie de cette contribution annule la deuxiéme 


partie de la contribution précédente, et il reste : 
la = 


Pp 
zs Sen (z, ©) )- 2 #(a2', 0") ula’, 2") IMer, x) Yale, 9)dedz’ 


hy 


Dea y SS 8b5*(x, 2) of | lc} L(x, o)dx 


eras 


-yJ x by.*(z, mee tN Sie acer 
$, t= 2. =1 


32 APPLICATION DE LA METHODE DU CHAMP SELF-CONSISTENT 


On en déduit pour jes autres types d’interaction : 


81g. = 


n 


Df Dae (x, s) y Yu*(x’, 9) -dy(x, o’)- J2(x, x') s(x’, «)dadz’ 


s,o/" s=1 Mo 1 





+f y v v Pats (x, ¢) Beall oe] dole odd 


$,t= Ma 
n 


Sj Sr *(z, c) 


bile Fg | du(a, 2) o)ax 





64143= 


xsz y ae 2) BY dala 2) dale 2) Ce, 2!) te, odode! 


+Ofdy Dd) eye*( (x, ¢) fel [et wate, eax 


s,t=1M=1 


be P 


kf ND averted fe] TSM] ele, oar 


o1y4= 
Y | Seve, 9. y D, Yu (a’, 0!) Yue, 0) Tylor, 2") tela’, Vdd’. 
g,0 ro uM. 1 


4 
: bs(x, o)dx 


np 
+f dS Dewey 
s,t=1M=1 


c 








“Y 


P 


-rf> SD eyez, 0) |! 


s,t=1M=1 





4 
vA V(x, o)dax. 
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Valeur de 0,1. 
Nous pouvons maintenant rassembler la valeur de 6,1. 


2) 41= > f Youre “|. = E)-t, 9) 
a” s=4 


+ [Ht +S Htelba(e, 2) — SY Hastale, 2) 


i k=1 





















































P ; Yu* (a, 0°) +$u(a', o')- Ja, 2')-Ys(a, 2) 
=e Vu*( 4 a’) + Yu(x, 3’). - J3( D x')-Y4(2', 3) 
sf py | +e 2’) sale ¢): 1%, x’) -Ys(x, 3’) ms 
ena bate o)-TMes 2) 4400) 
/ 2 3 J4iN 
eS ll at al Be 
atin! [— 3 a ae ee i Pate t |J8 ae M J Jae,2) | 
-dx. 
Valeur de 9,1. 
La variation 0,1 s’obtient facilement par la méme méthode 
That ale é¥,*[B, —E+L,+C+ D4] Yedr 
Sn+1,°**Sntp 
ey. 
~2 4 2 abt (a, +] (By — E) a(x, 2) + 
go) gt eee = y Hxx|Yi(2, 7) — y Hxuhx(2, 2) 
K=1 a4 
+ Cv,(zx, ¢) 
Ym*(2’, 0’) + Yn(2x’, a’) -J he fe u(x, fa 
Um*(x', 5’) + Ym(x, 9’)+ (e:% (2, 5 
+E sd! eee imi 3) J3(x, x) -Wy(z, “)( 
peso tt Latmles 2) 9 2) oy A T \J4im } 
fe P| Ii", oH ha +|" za Be ieed nine ml 4 Jee 
FP LIT bm) UH mI ate 
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Equations des fonctions d’onde partielles. 

Pour obtenir les équations des fonctions d’onde individuelles, 
il faut égaler A zéro le coefficient des variations indépendantes 
dans 0,1 et 3,1. Les variations totales 3,1 et 3,1 étant de la forme 


6,1 -> » bY,(x, o)[-- + ]da 
fey s=—1 


on voit que les équations des fonctions d’onde partielles s’ob- 
tiennent en annulant la quantité entre crochets. 
Il est inutile d’écrire ces équations sous leur forme actuelle, 


5. Valeur de I’énergie. — Considérons l’équation de la fonction 
L(x, a). Cette équation n’est pas écrite, mais il suffit de supposer 
égalée a zéro la quantité entre crochets de l’expression (33). 

Multiplions & gauche cette équation par vi(z, o) et intégrons 
dans le volume du noyau, en sommant sur les deux valeurs du 
spin. On obtient en utilisant les propriétés d’orthogonalité et 
de normalité des fonctions / une relation qui détermine E. 

P 
B,—E+ > H«n4+C 
K=1 


+ Dll al Pod + fol” 


" y y {T, rae P 
m=1T=1 


j2 








ie ool nee a aa 
mith ml lett ml |r] 
—S [Fall + Pt + ol hat + ll = 


L’énergie a donc pour valeur 

















j2 




















» P 
(34) E=) H+ > We 
k=1 K=1 


+ SD DU lel tel O14] 
Hele tal |e 
+] |" i Ss | 
+{%,|° FTA ral 


np 
SSP Pala ble Pal 








| 
i 
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r"| es ee 
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Les premiers termes représentent Pénergie cinétique. La 
deuxiéme série de termes représente énergie potentielle entre 
particules semblables (dans notre cas entre neutrons). Les termes 
en G! représentent l’énergie due A une force ordinaire, la partie 
positive correspondant a l’énergie habituelle et la partie négative 
a la correction d’échange due a Pantisymétrie de la fonction d’onde 
(14) vis-a-vis des particules semblables. Les termes en G4 repré- 
sentent l’énergie due a une force d’échange du type de Heisenberg ; 
on voit qu’ils ont la méme structure que les termes en G!, mais un 
signe opposé. Par conséquent l’action normale d’une force de Hei- 
senberg est identique a l’action d’échange d’une force ordinaire et 
inversement. La méme réciprocité existe entre les termes en G2 relatifs 
ad une action d’échange de Majorana et les termes en G' relatifs a 
une force avec échange du spin. 

On vérifie aisément que cette valeur de l’énergie satisfait a la 
relation importante 


(35) E => [ vena: 
=B,+B,+C +3 [vee 


En portant cette valeur de E dans les relations (32) et (33) on voit 
que les coefficients de s(x, c) et t.*(x, c) se détruisent, a part les 
Laplaciens H! et H? et les équations des fonctions d’onde indivi- 
duelles prennent la forme plus simple 


(36) H?-4,(x, 2) 
Ye*(x', 9')-4i(2', o’)-Gl(z, x’) -¥s(z, 5) 
— ¥*(2', 3’) )-4e(z, 3) Gl(x, x x). “Y,(2', 3’) dz’ 
+E y + HM a’, o')- Ya, 3’) G(x, x')-¥,(2', 2) 


tas 


Ch 0 a eed O28 6 le WO Oe 6 626 LG pe A109 70 Oe 6 08.68 


Ya (x', 3’) + bu(z’, 0’) J*(x, x')-Ye(z, 9) 
4% Yu*(z', a’). Yur (x, a’) J*(z, 2’) V(x", 3) 


+Ef i). 





= >> es'ba(Z, ¢) 


k=1 
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(37) H?.d4(2, 3) 
5 Ym*(x', 5!) Ym(x', 0) I*(ax, a) - Ya (2, 2) 
+d i Dy fF bnt(o’, 0") -Yn(a, 2) I(x, 2')-ta(a, 0) § de’ 
a! m=1 +b a\¢ 0p «4 w.eK6. 5.0 apie nie we ace ate. hie = ls aia aharees 
P 
= > Anu Yx(2, 9) 
K=1 
D’apreés 
n : L 
Wes py mf 





on voit que 


Il en est de méme pour ss. 
On peut effectuer sur les fonctions ) une substitution linéaire 
orthogonale de fagon a rendre diagonales les matrices Aix et ty. 
En effet, les n? coefficients qui définissent une substitution li- 
néaire sur un groupe de n fonctions vs par exemple, doivent 
satisfaire aux conditions suivantes : 
n(n — 1) 
2 
n(n — 1) 
2 


conditions Aw = O pour i x k 
conditions d’orthogonalité des nouvelles fone- 
tions ts 


n conditions de normalité des nouvelles fonctions ys, soit au 
total n? conditions qui déterminent les n? coefficients. 





6. Definition des densités mixtes et de divers opérateurs. Dé- 
sormais nous rétablirons la symétrie dans les résultats en intro- 
duisant les opérateurs G},.. G4 d’interaction entre protons. 

Nous posons par définition 


| (219 | en | Xg4) bm™(X1, 5) Yn(Xg, Tq) 


=. 
anol 
m=1 
P 
| (27492 | Pp | %y%4) = SI Yu*(2y, %)-Ya(g, 24) 
\ M=1 


D’autre part ,en désignant par P un opérateur de permutation 
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des coordonnées d’espace, par Q un opérateur de permutations 
des coordonnées de spin, nous définirons les opérateurs intégraux 


Gri, GP, G,FQ, n'PQ 

Gri PQ, Gr?Q, G,*P, Grit et les termes 
J/1, td 5,°0; J,*PQ analogues en 
Jn1, tS n°Q, w@PQ G,1:--Gp* 


par des relations telles que 


Gut -4(c,2) =D,  YmPle' =!) Galle, 2) Yale e) Hee) 
g! M=1 


— 2 
G4PQ-u(c, 2) =D) [DS ete Gla, 2) tule, 9) He od 


o “ M=1 


P 
,°Q-4(2, 2) = > if S. gate 2") Je, 2!) ula’) Ya, ode 
o! M=1 


Jn*PQ: Y(x, 3) >; >? f » Ym*(2', 3) : J*(x, x’) : Ym(2, 3) : f(z’, a’) de 
o m=1 


On peut également introduire les densités (38) ce qui donne les 


définitions équivalentes 


Pee Star PotD, | o-Oe< Ire ee) Hite 
aver: Qla’a! |e a’s')-¥(x, 0) = as" |e | 2’) 412.0) 


7. Equations des fonctions d’onde partielles sous forme d’opéra- 
teur. — Avec ces définitions les équations (36) et (37) prennent 


la forme simple : 
(40) | Ht + Gal(1 — PQ) + Ga(P —Q) + GaX(Q— P) + Ga(PQ — 1) 
4 J + IpP + JQ + IptPQ{ Hila, 2) = M(e,9) 
(41) \H? + G,'(41 — PQ) + GA(P — Q) + G,8(Q— P) + Gp(PQ— 1) 
+ J,1 + JP + JQ + Jn*PQ | bu(v,2) = uba(Z, 2). 

Les opérateurs entre accolades dépendent des fonctions v et 
par conséquent le probléme n’est pas linéaire. Mais les notations 
(39) montrent que cette dépendance des opérateuss s’effectue par 
Vintermédiaire des densités p. 

Par conséquent dans (40) opérateur entre accolades est unique. 
Quel que soit indice de la fonction 4 considérée ct 11 en est de 
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méme pour (41). On voit done que les fonctions 4 sont toutes 
fonctions propres d’un méme opérateur et la supposition que nous 
avons faite initialement sur l’orthogonalité de ces fonctions 
n’a introduit aucune contradiction dans la théorie. 

De méme, les fonctions Yy sont un ensemble de p fonctions 
prises parmi les fonctions propres, en général en nombre infini, 
de l’opérateur défini en (41). 

Les opérateurs, P, Q,... etc., qui ont été définis en (39) comme 
agissant sur une densité p et sur une fonction / peuvent aussi 
agir sur deux densités de la facon suivante 


P(x | pp | xs)-(x'o" | gn | 2’o’) = (x0 | pp | 2's)(2'o’ | on | xo’) ete... 


8. Autre expression de Il’énergie. — Avec ces opérateurs et les 
densités p ,l’énergie E peut se mettre sous la forme plus explicite : 


n Pp 
E= a+ >} ay 
k=1 K = 4 





vf 
ss 3D / [Gat(x,2')(1 — PQ) + G,2(z,2’)(P — Q) 


+ Gri(x,2')(Q — P) + Gaf(x, 2')(PQ — 1)] 
(2 | pm | ae)-(a,0" | pn | a'o')dor dar’ 


(42) +f [J(a,2’) + J%x,2')P + J%(a,2')Q + Ie,2')PQ] 


(x2 | Pn | e)-(a's' | op | a's’ )dada’. 
1 
+ Daf [Gp'(x, x')(1 — PQ) + G,*(x, 2')(P — Q) 


+ GpX(x, x')(Q — P) + Gp(a,a’)(PQ — 4)] 
(12 | pp | o)-(2'e! | pp | 2's" )dx da’. 
= Eein.né + Ecin.p +I Enn =e Enp =e Epp 


9. Correction due au mouvement du centre de gravité (#1), (84), 
— Les axes auxquels sont rapportées les fonctions d’onde étant 
quelconques il peut y avoir un mouvement du centre de gravité, 
et comme nous cherchons l’énergie de liaison du noyau, il convient 
de retrancher de la valeur (34) de l’énergie. l’énergie cinétique du 
centre de gravité. 


Si on désigne par P la quantité de mouvement du centre de 
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gravité et par M la masse du noyau, énergie cinétique a pour 


valeur en mécanique classique 


1 
Da 
a P°: 


A cette expression correspond un opérateur qui permet de cal- 
culer la valeur moyenne de l’énergie cinétique du centre de gra- 


vité : 


1 
Eg => fv om Powe 


is] 


avec : 
2 


i=1 
n+p n+p 
=> pet > per 
Sins ifm 
ie tj 


Les opérateurs pi et pj ont la signification habituelle, ce qui 


donne 
n+p n+p : 
2 2? 3? a? a? P) ) 
(43) P? -|> (= - aye + =) a > (— + Syays = sme) 
L=4 ane 
h2 n+p n+p 
K—-pa= Dat D Gu 
i=1 ij=4 
iF] 
L’opérateur Gi = 2m,H} = 2m,H? donne naissance aux 
termes 
m Pp 
mM Ms | 
(44) M ») Hi + Wy mai 
= 


i=1 
Si on néglige la différence entre la masse du neutron et celle du 


proton 
m= m,= mM 
M=Nm 
N=n+p 


cette correction a pour valeur 


1 
N Kein. 
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D’autre part, si on pose 


Gij 
2M 





— Ky 


cet opérateur agissant sur deux particules est analogue aux opé- 
rateurs G,(ai, 2), J,(ai, x), et on obtient les termes 


(45) > > (Hl) fle] 
a3 SIA 


+S Sal a}—l [ef] 


Pour expliciter ces formules il faudra se souvenir que ]’opéra- 
teur K a en définitive la valeur suivante : 


m| 


h2 22 32 a2 
oor ~ 8x2M (a 3 dy: oy’ + aa?) 


Nous avons vu que l’énergie E satisfait a la relation (35) 


eS [ wtH Wd: 
Gyros 


Pour que cette relation subsiste quand on effectue la correction 
pour le mouvement du centre de gravité, il suffit de retrancher 


Popératenr 5, P? de l’Hamiltonien (1) 


Or: 
i) ie Ms a 
ene aoa | YY 
P= yf Het OS ae+ » Ky. 
i—4 i=4 i,j=1 
tAy 


Le dernier terme s’écrit également : 


nm n+p n+p 


2 y ky t2) VY xy+a2 DHEA, 


i>j=1 i=1j=n41 i>janqs 
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L’Hamiltonien du probléme devient donc : 


Snes 
(46) H=— gm, (1M) Sate (1B) S A; 
= i=n+1 


ord 34 , 
mS > | G,(r ix) — 2K, + G,(r iz) Pie + G3(r. in Pie + Gi(r in) Pie | 


t>k=1 
n+p 
Te > [ Go(re)—2K e+ G; (rie) Pie +Go(r *\ ps r+ Grin) Pig | 
i<k=n+1 
n n+p 
+>) > [Pra)—2Kat-F(ru)PR+ I (ra)P2+ I(r) PE]: 
i=1k=n+1 


Si on était parti de cet Hamiltonien (46) on aurait trouvé la 
valeur de l’énergie corrigée du mouvement du centre de gravité, 
mais les équations (40) et (41) auraient été compliquées. Il faudrait 
remplacer 


H?! par (1 — it) H! 


H? par (1 — ") ue 
Gr) par Gp — 2Kn 


Jn par Jn — 2Kx 
Jp! par Jp' — 2K. 


Les opérateurs Kn et Kp ayant les significations suivantes 


Kn}(z, 3) Le x Ym*(2', 3')K(z, x’) + Ym(2' , a’) +(x, s)dx' 
Kp: (2, I=L Dore K(a, 2')«du(x’, 0’) U(x, 9)dax 

ty af hm BY , Om dy 
--in D Lf (2', y, 2,2): | ee ee 


dx'dy' dz’. 
On voit que les équations que nous allons obtenir, quoique for- 
mellement aussi simples que (40) et (41), sont en réalité beaucoup 


plus compliquées car l’opérateur différentiel ne se réduit pas a 
un simple Laplacien comme H?! ou H?, mais il comprend de plus 
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des dérivées du premier ordre figurant dans les opérateurs inte- 
graux de l’équation. Pour cette raison nous conserverons les 
équations (40) et (41) comme les équations du champ self-con- 
sistent, tout en considérant l’expression (46) comme étant |’ Hamil- 
tonien réel du probleme. 

Par conséquent les fonctions {i et yy déterminées par (40) et 
(41) ne sont pas celles qui introduites dans (15) et (16) donnent la 
solution la plus approchée. Pour obtenir des fonctions meilleures que 
di et Yy il faudrait intégrer des équations telles que (40) et (41) 
mais corrigées du mouvement du centre de gravité. Nous venons 
de voir que ces équations modifiées sont beaucoup plus compli- 
quées que les équations primitives. Aussi nous contenterons-nous 
d’un moyen terme en faisant dans (40) et (41) la seule correction 
qui n’altére pas le caractére de ces équations : 


le remplacement de H? par(4 — m2) H? et celui de H? par 


PS Nh 
( mt ) He, 

Avec lHamiltonien (46) l’énergie E donnée par (42) et corrigée 
de (44) et (45) satisfera toujours a la relation (35). 


10. Perfectionnement des résultats au moyen d’un calcul de 
perturbations (!”). — Suivant la méthode de Moller et Plesset (27) 
nous considérerons la fonction d’onde approximative (14), (45), 
(16) comme étant la solution exacte d’un probleme dont nous 
désignerons l’hamiltonien par %. Nous appliquerons a cet hamil- 
tonien une perturbation ayant pour valeur la différence entre 
Phamiltonien exact H de notre probléme nucléaire défini en (46) 
et ’hamiltonien 4 de la solution approchée (14). 

L’hamiltonien # a pour valeur une expression de la forme 


n n+p 
berate S #?+C 
i=1 t=n+1 


dans laquelle 46! n’est autre que Popérateur de gauche de Péqua- 
tion (40) et 62 celui de Péquation (41). Les indices signifient 
que ces opérateurs agissent seulement sur les fonctions des 
variables 24, ci, de la jéme particule. Dans tous ces calculs nous 
considérons un état bien défini du noyau, caractérisé par les n 
fonctions et les p fonctions x ainsi que par la valeur E de 
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Pénergie qu’on en déduit au moyen de (42), (59 bis) et (60). 
Cet état sera l'état fondamental si les n fonctions 4: d’une part 
et les p fonctions vy d’autre part sont respectivement les n pre- 
miéres fonctions propres de l’opérateur 44} et les p premiéres 
fonctions propres de l’opérateur #6?. Ce sera un état excité dans 
le cas contraire, Mais dans tous les cas le calcul de perturbation a 
pour but de perfectionner les résultats relatifs A un état bien déter- 
miné du noyau caractérisé par les n fonctions i les p fonctions 
vy ainsi que par l’énergie E et les opérateurs 46} et 46? formés: 
avec ces fonctions. C’est une constante ajustée de telle sorte que 
la premiére valeur propre de l’opérateur 4 soit égale al’énergie E 
définie en (42) et diminuée de (44) et (45). Or cette premiére valeur 
propre a pour expression 


n n+p 
= > {i > a kae 3 v0 | Wdt + C. 
Ty" Onip i=1 i=n+1 


Comme l’opérateur entre crochets se décompose en une somme 
d’opérateurs n’agissant 4 la fois que sur une sewie particule, on 
trouve facilement, en répétant le raisonnement qui nous a con- 
duit aux expressions (21), que 


P 
4 ->/ Fuster) Hn Y;(x,9) dx + Quits )-Hn*b_(x,5) >) dx + C 
¢ ya | K=1 
n P 
= Hxu+ Yu +C. 
i=1 K=1 


Comme nous voulons que 
n=E 


on trouve 


n P 
C=E—xu— » Ax. 
i=1 K=1 


Nous appliquerons la méthode de Schrédinger car nous suppo- 
sons que la perturbation est petite. Nous supposerons quil n’y a 
pas de dégénérescence. 

Dans ces conditions & chaque solution approchée ’ correspond 
une solution exacte ® qui en est voisine, et il en est de méme 
pour les valeurs propres approchées E et exactes «. 
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On peut done développer ® et ¢ en série de la facon suivante 
(b= WV +6,4286,4... 
( e= E+ de, + ie, 4. 
En introduisant ces développements dans Péquation 
(H%+2£)b=ceb avee P=—H—HK, 
et en identifiant les coefficients des m4mes puissances de i on ob- 
tient : 
(48) HV = EW 
(49) (#6 — BE), = (e, — ayy. 
(50) (4% — E)®, =e, + (e, — £)4,. 

L’équation (48) est identiquement satisfaite. 

On démontre que l’équation (49) n’a de solutions non nulles 
que sile membre de droite est orthogonal aux fonctions propres 
de Popérateur de gauche égalé a zéro. Or Popérateur J6-E n’a pas 
d’autre fonction propre que ¥, donc 





y [vt —swae=0 


T° . *on+p 


(51) g= a wed: = 0, 


car, d’une part 
€=H—H# 


Df verwa: = E, 


et d’autre part nous avons choisj Vhamiltonien % de telle sorte 
que sa premiére valeur propre soit 


DY | vse: = E. 
Tyeee 


Ainsi, la correction du premier ordre de énergie est nulle, mais 
cela ne prouve nullement que l’approximation réalisée soit bonne : 
cela est di uniquement A la fagon dont nous avons choisi H#. 
On obtient de méme a partir de (50) et en utilisant (51) : 


Ey => DY f vevnye. 
Tyree 


et 
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Considérons maintenant la suite ¥’,, V’,,... des fonctions propres 
de Vopérateur J. 


Elles sont reliées aux valeurs propres correspondantes par les 
relations 
(52) HOW; = Ei. 

Les valeurs de l’indice z sont choisies de telle sorte que la valeur 
zéro corresponde a la premiére approximation ’ de la fonction 
d’onde, et a la valeur E de l’énergie considérée jusqu’a présent. 

Les fonctions ¥; formant par hypothése une suite compléte de 


fonctions orthogonales, on peut développer ®, et ®, suivant ces 
fonctions 


(53) b, = Vary; 
i 
>, = ~~ a?zW;. 


En portant (53) dans (49) et en tenant compte de (52) on trouve 





s¥ a}(E; — E,)¥i = — £%o, 
i 
dou 
>2 Jf voor: 
54 ws ayers . 
(54) Ga = —7 7 pour i + 0. 
ag: = 0. 


En tenant compte des expressions (53) et (54) la correction ¢, 
de l’énergie au deuxiéme ordre s’écrit : 


ee Degg d frome E [rsenae 
i=+0 atte ol jeee 
Pour expliciter les calculs il faut définir comment nous for- 
mons la suite des fonctions Vi a partir des fonctions propres des 
opérateurs 4, et #,. Remarquons, (*’), qu’en plus des n fonctions 
va qui entrent dans le déterminant ¥’,, il existe une suite infinie 
de fonctions propres ¢2 de l’opérateur 4d et il en est de méme pour 
les p fonctions Ye qui entrent dans le determinant ’, et aux- 
quelles il faut adjoindre une suite illimitée or de fonctions propres 
de Yopérateur 4’. 
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Les opérateurs d6' et J6? sont construits un fois pour toute avec 
les fonctions {i et Yu et les fonctions 9 sont précisément celles de 
leurs fonctions propres qui n’entrent pas dans leur définition. 

Pour obtenir toutes les fonctions W, il suffira de remplacer dans 
(15) certaines fonctions ta par des fonctions 9,, et de méme dans 
(16) certaines fonctions 4 par des fonctions Op. 

Une telle fonction sera désignée par 


Les indices a peuvent prendre toutes les valeurs de 1 4 n, les 
indices c de 4 a p, et les indices « et I, de 4 a V’infini. 

Avec ces notations les éléments de matrice qui entrent dans 
(54) et (55) s’écrivent 


CD ;TA 
» furs 14g £V ode 
a 


car lopérateur # se décomposant en une somme de termes qui 
agissent au plus sur deux variables a la fois, toutes les fonctions 
¥i qui se déduisent de YW, par l’échange de plus de deux fonctions 
 donnent une contribution nulle a cause de Porthogonalité des 
fonctions ¢. 

En définitive la correction du second ordre a énergie s’écrit : 


Ee ! CD ; ra 
aps Greens) SpfieePoe ees 


ba CD i eae 
aBlA er eye ay 
GD TAC = 
oe fveree OP 5 ds 


en posant : 
c 
Ap =D f ertaetepa 
c 


et en tenant compte de (47). 


11. Décomposition de chaque onde individuelle en un produit 
d’une onde d’espace par une onde de spin (25), (7°) he spin, ne 
. . . : 
figurant pas explicitement dans Phamiltomen, nous pouvons sé. 
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parer dans chaque fonction d’onde individuelle, une fonction de 
spin et une fonction d’espace 

Yala, 9) = 94(2) X 24(9) 
Le nombre quantique de spin x n’est susceptible de prendre que 


2 
correspondantes par les notations 


a(a) eb (2). 
+ —— 


1 1 i ; 
deux valeurs + 5 et —— et nous désignerons les fonctions d’onde 


La variable de spin o ne peut prendre aussi que les deux valeurs 


1 
+5. 
ppl. 

Les conditions d’orthogonalité et de normalisation se scindent 
en deux groupes de relations 


if gi* (x) -ex(x)da = Cie 
SY ai*(2)-24(2) = biee 


¢ 
Nous supposerons tout d’abord les fonctions i rangées dans un 


ordre tel que pour : 


O<cicn a4(s) = (9) 
m<aign 3 ea, 
0<K<p ney) 
pa<K<pP 0x(3) = 2(s) 


_— 


Cherchons comment s’exprime, avec Ces nouvelles hypothéses, 


énergie (34). 
Les termes d’énergie cinétique ne changent pas de forme, car : 


Ha = ~ df e*(x, 0) H?-Ya(z, o)dx 
— Si xe*(°) - ap () . fe) . H1e,(x)dx 


et comme nous allons faire sortir la variable de spin de toutes nos 
écritures nous pouvons transposer toutes les notations introduites 
jusqu’a present et les appliquer 4 des fonctions d’onde ne conte- 
nant que les coordonnées d’espace. 
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Nous posons donc 


Hi, = | ge*(x) -H¥94(x)de 


(tk | G* | Im) = if ve*(x) + o4*(x')-G*(x,a')gx(2) Pm x')da- da. 


he 

S) ci*(2)- ex’) = (@ | em | 2’) 
i=1 

(56) 
y Fiaal ( %)+ g(x’) = (x | png | x’) avec Ny =n—n, 

i=n,+1 

Py 

D> %*(2)-9x(x’) = (2 | pps | 2’) 
K= 


P 
Di ex*(2)-9x(a') = (| pala’) avec: pp=p—p, 

\K=p,+1 
Calculons maintenant par exemple un terme d’énergie potentielle 





oe 








m=1T=4 
ee tig ce terme 


y x. a3 > Yn (2,6) + Y(t! 0") - J3(a2,a0')n( 22,0") « Yon ae" c)dxde' 


ar Cy eae 


Calculons a part les sommations sur les fonctions de spin 


ee i 
XY Lex*(e)-an(s)- Yan*(o') an 


m=1T=1 ¢ o! 
- EP Ero. Vege 
m=1T=1 o¢ 


= i > y Lerle)-a( (2). S) a*(e').a(0’) 
m=1T= Be 3 gat, = 


+ S; X de a* (a (2). §} -a*(q') -a(o' 


m= net Tt a a! 


++ y ) y Sy 2%(2) -2() ‘ 83 -a¥(a')-9(3') 


m=n,+1T=p,+1 ¢ o! 
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A cause de l’orthogonalité des fonctions « les seules contribu- 
tions non nulles sont celles relatives aux sommations 


nt Py n Pp 
pi NED) 
m=1 T=1 m=n+1 T=p,+1 
On obtient donc le résultat 


( » y a >» y ) Pot (ar) « a(x) - Pm*(2') - Pm(2') 


m=1T=1 m=n,+1T=p,-+1 
J3(x,2') dada’ 


= { ley l2-( 2 | Pn, | 2"): S(a,2') dada’ 
aE f (elep, 2): (2' | Pn, | 2’): ): J3(x,a") «dar de’. 


Les autres termes se calculent aussi facilement et on obtient 
expression définitive de Pénergie 


w+ DW 


(a | on | x)+(x’ | pn | 2") -[Gu'(a,2") — Gu¥(a,a') |drdar’ 


(57) 


fi 


Dole it 3 


+5 
x7 1 (x | em | 2!) | 2-[Gn2(x,2") — GnP(x,2') Jd der 

+f len 120-C'Len,|2Goe2) — Gita) d 
24 5 ff (x | Pn, |2)-(2" | Pn, | 2 Ga?(a,2") — Gu2(ar a’) \dacdar' 
ot : fi (| en, |’) [*LGn*(2,2') — Gal(a,x’) |dx- da’ 

+5 f Ue beng 2) Late) — Cae ede 

Is 4} (x | pp | 2) +(x! | pn | 2’): 3*(x,2')dvde 

+ Ak (x | pp | a")-(a! | pn | a) - J*(«,2")dardar’ 

+ fc (| ep, |2)-(2" Len, | 2") S(e,a') - deed’ 


+ f (el epyl2)-(2' lng | 2) Pee) ede 


4 
M. MarTricon 
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(57, suite) +f ( | Pp, | a’) «(2 | Pn, | x) - J4(a,2')dada' 


+ fe lepyL2)-(e! Leng |2)- IM! ded 
+ termes d’énergie potentielle correspondant A G,!,---Gp*. 
Pour obtenir l’énergie de liaison du noyau nous devons retran- 
cher de l’énergie totale (57), énergie cinétique moyenne du centre 
de gravité. 
Avec les nouvelles hypothéses sur les fonctions de spin, cette 
énergie moyenne a pour valeur : 


" P 
Mia 1 Mey Ae 
we a Het + yp 

i=1 i=1 
( my ny Ny+Ng Ny+N, 


yys+ B+Y Jr) — ae 4 


t=1k=1 i=n,+1k=n,4+1 
Ny+Nny 
423) SN (ik| Kl ik) 
i=1k=n,+1 
n P 
+2 HM (Tm| K| Tm) 
m=1T=1 


Piper PitP2 PitP2 
+($S4 er i )t [((TM | K | TM) + (TM | K | MT)] 


M=1T=1 M=p,+1T=p,+14 


Py PitP2 
+2) §} (™|K|TM). 
M=1T=p,+1 


L’opérateur K a pour expression 
i a F) o 
me sani se-ay 15 ay oy + saa) 


Nous remarquerons que cet opérateur agit toujours sur le pro- 


duit d’une fonction de z, y, z, par une fonction de 2’ hs) pie 
Ore 


onc Tal Na Ort eee 


On peut donc le mettre sous la forme 


h? PY a PY 3 a Py 
K = — ars “ae ag tay “ay +E XB) 
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Nous sayons que tous les noyaux connus a partir de l Hélium 
satisfont a la relation 


pan 
Pour permettre le développement des calculs nous ajouterons : 
Pi P2 Sm <M 


On peut faire d’autres hypothéses sur la grandeur relative des 
nombres Pj, Ps; M4, M2, mais cela ne change rien a la fagon de con- 
duire les calculs. 

Comme le principe de Pauli nous permet d’attribuer la méme 
fonction d’onde spatiale a deux particules de méme nature mais 
de spin différent, nous poserons 


¢4(2) = Pi4 n,(*)s 


et les densités mixtes auront pour expression 


SY) ee* (x) -gelz’) = (| Pn, 1) 
(58) gh 
| SY} ex* (2) ex!) = (| Pp, 12): 


Ge 


D’autre part comme les fonctions de spin sont supposees données, 
les variations ne doivent porter que sur les fonctions d’onde spa- 
tiales, et par conséquent 


BYi(x, 3) = aa(s) X ¢a()- 


Par conséquent lorsqu’on considere le coefficient d’une variation 
pour chercher une équation d’onde, il faut comprendre dans ce 
coefficient la fonction «(z) et effectuer la sommation sur les deux 


valeurs de o. 
Considérons par exemple une valeur de l’indice i comprise entre 
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0 et n,: la fonction de spin correspondante 4a la valeur 2 (a), et 
le coefficient de la variation d9,(x) s’écrit, d’aprés (36) 


H19;(z) 


Fa fe | en |") -[Gu2(x, 2”) — Gaba, 2')}ea( ard 

+ [|p| 2)-{Gi%a, 2") — CaM, 2p) ae 

ke dp (2' | Pn, | 2')-[Gu(a, x) — Gar(ax, 2’) }ei(a) de’ 
me [ (2' | py, | 2) -[Gn8(x, ar’) — Gal(ar, 2')] + 4(22\der' 
+ J “Ca |¢p | 2')-32, 2!) -ge(a)de 

+ J bop |2) Je, 2) -94e!yde’ 

+ [bp L2)-I%e, 2) ula) 


+f (Le, La) Ia, 2) ade 
— 3 Aeiz a(x). 
k=1 


+ termes relatifs aux protons. 


Or, 4 toute valeur de i comprise entre 0 et n, correspond deux 
fois la méme fonction d’onde: 


g(t) eb ign (t) = ala). 


I] s’ensuit que l’on ne doit pas annuler séparément les coeffi- 
cients de 


seta) et byt, (x), 


puisque ces variations, étant identiques, ne sont pas indépendantes. 
On doit annuler la somme de leurs coefficients. 

On peut pousser la simplification plus loin et admettre que les 
fonctions d’onde spatiales des protons et des neutrons sont les 
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mémes. Avec cette hypothése les équations prennent la forme : 
(59) (aiH* + biH?)¢i(x) 


‘ gx*(a’)-92(x!) (Gut — Gat) gil) ) ,, 
a “( ~ +x s + ¢e* (a) x(x) = (Ga® — Gn®) -ge(2’) ae 
2 \( gala!) -gele)-(Gut — Gut) eile’) ) 
f(ad+ud eee *(x!)-gx(a!) (Gu? — Ga’) -¢ watt 
+ fi +3 )+0(S+) ‘) 
k=1 k=1 1 


+f (mH +i Send +n¥ ) 


band Just kent 
x')+ox(x')- I(x, is ie hae 
x')+ex(x)- J4(a, 2’) x’) 


eat(a!)-eu(2) (Gp — Gy4)-e!) 5 
“eked Ve H(2')-(Gy? — G,) gala) 


ooo oo PF OrPorPrPPrPo YS 


2 
2 
2 
1 
1 
2 
2 
1 
1 
1 
1 
2 
1 
1 


ORR ORR REP NP EP NPD 
ececoocorFRONKF KNOWN 


—) 
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On peut classer les opérateurs intégraux qui figurent dans ces 
équations en deux catégories : ceux dans lesquels la fonction 9, 
portant sur la variable « reste sous le signe d’intégration, et ceux 
dans lesquels la fonction gi portant sur la variable : x sort du signe 
d’intégration. 

Les premiers opérateurs peuvent se mettre sous la forme abrégée 


DY Qie(z)-¢2(2). 
k 


Nous excepterons de cette somme le terme correspondant ak = 1, 
et nous le joindrons a la deuxiéme catégorie d’intégrales. Nous 
désignerons l’ensemble de cette deuxiéme catégorie par 

Qui(x) - p:(). 

Pour rappeler que dans les sommations sur les valeurs de k la 
valeur k = i est exceptée, nous accentuerons le signe 2. Nous 
ferons de méme pour les Ai. 

Posons encore ; 
aiH} + 5;H? = cA. 
les équations des fonctions d’onde ont donc la forme 


(60) ciAgi(x )+ Qi(x x): - ¢4(2) — heipi(x a Mpe(Z) ) —S Qala) -ga(a). 
k 


12. Equations des fonctions d’onde satisfaisant a la répartition 
angulaire en fonctions sphériques (*) et correction pour I’énergie 
cinétique du centre de gravité dans ce cas. — Considérons par 
exemple un noyau caractérisé par 5 fonctions d’ondes spatiales. 

La variation que nous devons annuler a la forme 


(61) J Pedra oe HE bp4*(x)Agdx 
+f erst (x)-Agdax +f toate -A,dx 


Dans cetve expression As représente le premier membre de l’équa- 
tion qui se déduit de (60) en faisant tout passer dans ce premier 


membre, soit : 
Ag = GiAgi(x) + Qui(x) pil) — raepi(x) 


= hivpe(ae) + SY a x) 9x(2). 
k 


Nous supposons que chaque fonction (x) = g(r, 4, g) est le 
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produit d’une fonction inconnue du rayon par une fonction sphé- 
rique Y(6, 4). 
(62) vl, 9, %) = gna(r)-Yim(9, ¢) = gy(r)+¥p(9, ¢) 
en désignant selon la notation habituelle par n le nombre quantique 
total, par / le nombre quantique radial et par m le nombre quan- 
tique magnétique. 

La correspondance entre les valeurs des indices est indiquée 
dans le tableau ci-dessous : 





Les cing premiéres fonctions ont donc la forme: 
re 
y= £1(r) x Y,(6 o) = g,(r)- \/ ez 


Po = 82(I") X ro 


3 
eee p 


: 3 
24 = g(r) - sin 0-e'? 3- 


pee eae 
?5 = g(r): sin 9-e ? Be 


En portant ces fonctions dans l’expression (61) et en égalant a 


zéro on obtient : 


cy 
J eeittey-rtr [ oA, sin 0ediede-y / Ze 
Re 
me [eeyr(r)-r2ar fA, sin §d0-do- \/ iz 
| s 
| Az cos 0-4 / 7 | 
2 A, sin 6 —? 3 | sin 0-dd-do = 0. 
+ [testi dr | «+ Aqsin 6-e ; an 
| in 9-e'? 3 | 
\ + A; sin -é \/ & 
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En annulant les coefficients des trois variations indépendantes 
on obtient les équations qui déterminent g,(r), g,(r) et g,(7). 
Explicitons par exemple l’équation de g,(r) 


fa sin 6-d0-do-Y,*(0, 9) = 0 
63) i [Cd + Quilt, ® 2)-u(r)-Ya(0, ¢) 
—ar-gulr)Ya(0, 9) a daage(r)“Yu(9, ¢) 
+ YY) Quer, 6, ¢)-gelr)-Yu(0, 2)1Ys*(0, ¢) sin odode = 0. 
k 


Nous avons conservé la fonction Y,*(9, 9) bien qu’elle soit cons- 
tante pour indiquer la forme générale des équations. D’autre part 
Pindice k qui figure dans les sommations est purement symbolique 
et les véritables indices des fonctions g;(r) et Y.(9, ¢) correspon- 
dant a une valeur de & sont donnés dans le tableau de la page 53. 

Nous savons qu’en coordonnées polaires et avec Vhypothése 
(62) le laplacien a la forme: 


ce) O8i Yi 
Ao = Ee (7 se ea nee 1)g:| Y, 


r2 


Nous pouvons donc développer ainsi Péquation (63) 


(64) ¢, = ; 2 (n=) re ee s,| [yy sin 0-d0-de 
+ g(r) f Qulr, 9, ¢):¥,-Y,* sin 6-d0-de 
— g(r) [ Y,-Y,* sin 0-d0-d2 
—> hugulr) | Ye-Ya* sin 6-d)-d» 
= 


+) gx(r) [out 0, ©) Yx- Yi* sin )-d0-dg = 0 
k a 


En tenant compte des propriétés d’orthogonalité et de norma- 
lisation des fonctions sphériques, et de l’identité suivante, 


1a ag a2 
ey ads 2) a ee 
ror (- 7s) ~~ er? (rg) 
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on voit que l’équation (64) s’écrit : 


: Ul + 1 
o,[=4— iba fi| + Qulr)fi — sh 


=a: darfe 7>¥ Qu(r) fe = 0 
k k 


avec: 
fi(r) = r-gil(r) 


f A= An 


Qu(r) = | Qulr, 0, @)-Yy-¥y* sin 0-d0-dy 


Qual?) =| Qu(r, 9, ¢)* Ye: Yi* sin 6-d0-d¢, 
\ 


Avec ces nouvelles fonctions fi la condition d’orthogonalité et 
de normalisation des fonctions 9 devient : 


+ +2 
if gi*(r)-ge(r)-r2dr = i fi*(r)-fi(r)dr = Six. 
0 0 


Correction pour l’énergie cinétique du centre de gravité dans 
Vhypothése (62). — L’opérateur K se calcule facilement en coor- 
données polaires. 


87°M 





—ar Kk= 
: a2 . cos8-cosg 9a sin 9 2) we , 78 
6- .— + + .— — ——- — }-(sin® -cos¢ -s4 °°: 
(sin ome? art r 00 = =6r-sin§ o¢ ( ¥ ort ) 
: ; 3 cos§-sine 2 cos 9d oe ee SP 
° po — ————— -— es (AITO + 81N Or #5 
+ (sin 6 sin ¢ ate = a er =| ( 9 aoe ) 


Pee sin) a ae sin 0’ >) 
BM Or apc ap 30) No. Or’ ra’ 


Nous allons calculer un des éléments de matrice 
(ik |K| ik) qui se présentent dans la correction de 
l’énergie, en supposant qu’une fonction : gi par exemple est dans 
état s. Cette fonction ¢« ne dépend donc que de r, et Vélément de 
matrice s’écrit, 2 un facteur constant pres : 
f(r) -sin 9-008 ¢- SP resin 0-d)-ay X [ ceKen d? 


a 


Stod fa : > 
+ f get) -sinOosin SP rar. sin 0-d)-d x | oe* Ka! + Gear 


0° A PS r > 
+ [ sep) -00s0- FF tar sin 0-a-de x | on* Ka ox: dr. 
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en désignant par K’,, K’, et K’, respectivement la deuxiéme, la 
quatriéme et la sixiéme parenthéses qui figurent dans K. 

Les intégrales portant sur des variables non accentuées se décom- 
posent chacune en un produit d’une intégrale portant sur 7 par 
une intégrale portant sur 6 et 9, et les trois intégrales portant sur 
6 et » sont nulles, par conséquent 

(ik |K| tk) = 0 quand qm est dans l'état s. 

Pour que les termes d’échange (ik |K| ik) soient nuls, il faut 
que les deux fonctions gi et o soient dans état s. 

Dans le cas du noyau d’hélium qui ne contient que des parti- 
cules dans l’état s, la correction pour l’énergie cinétique du centre 
de gravité se réduit aux termes diagonaux 


“A 4: iy 7 Pi Pa 
3 ( > Ha + a ne + a ys Hiv? + y a) 
1 


i=1 i=1 i=1 


II 


INTEGRATION DES EQUATIONS DU CHAMP SELF-CON- 


SISTENT ET APPLICATION NUMERIQUE AU NOYAU 
D’ HELIUM. 


1. Généralités. — I] est A peu prés impossible d’intégrer ana- 
lytiquement d’une facon rigoureuse les équations intégro-diffé- 
rentielles du champ self-consistent. De méme que dans les pro- 
blémes atomiques, l’intégration rigoureuse ne peut se faire que par 
les procédés du calcul numérique. Aussi se contente-t-on la plupart 
du temps de solutions approchées. Pour utiliser cette méthode 
d’une facon approchée on utilise la méthode de Ritz décrite page 42. 
On se donne une suite de fonctions orthogonales satisfaisant aux 
conditions de continuité et aux conditions aux limites requises, 
et dépendant d’un ou plusieurs paramétres arbitraires ; avec Ces 
fonctions on forme les expressions (56) et (58) des densités et (57) 
de énergie ; on obtient ainsi une expression de Vénergie dépen- 
dant {d’un ou plusieur paramétres a, 6, etc., eb on détermine la 
valeur de ces paramétres par la condition de rendre cette énergie 
minimum, soit : 


aE ok 
Serb Rar Se 

Cette méthode a été appliquée primitivement par Heisen- 
berg (#6). Cet auteur a supposé simplement l’existence de forces 
d’échange de Majorana entre protons et neutrons et de forces ré- 
pulsives de Coulomb entre protons ; il a utilisé les fonctions 
propres de loscillateur harmonique, le paramétre 4 ajuster étant 
la fréquence de l’oscillateur. Tous les auteurs qui ont travaillé la 
méme question aprés lui ont adopteé la méme suite de fonctions 
orthogonales, mais en utilisant des lois de forces plus complexes, 
ou des fréquences différentes pour les fonctions d’onde des neu- 
trons, et les fonctions d’onde des protons. 
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Les résultats obtenus par cette méthode ne sont pas quantita- 
tivement bons : les énergies sont notablement trop élevées. Cer- 
tains auteurs (28) suppriment ce désaccord en ajoutant arbi- 
trairement A l’énergie calculée par les formules données précédem- 
ment, une fonction linéaire du nombre des particules dont ils 
ajustent les deux corfficients de fagon a retrouver exactement les 
énergies de liaison de deux noyaux déterminés. Mais cette fagon 
de procéder n’est pas justifiable théoriquement. Il est évident 
qu’en partie ces résultats proviennent de ce qu’une fonction 
d’onde en forme de déterminant ne peut représenter correctement 
Pétat d’un systéme dans lequel les corpuscules voisins sont en 
interaction trés profonde les uns avec les autres. 

I] est quand méme intéressant de savoir quelle est la part du 
désaccord qui est attribuable au principe méme de la méthode du 
champ self-consistent, et celle qui est attribuable au fait qu’on n’a 
appliqué cette méthode que d’une maniére approchée. 

Nous nous proposons donc d’intégrer numériquement |’équa- 
tion du noyau d’hélium et de comparer la valeur de l’énergie ainsi 
obtenue avec celle que donne la méthode de Ritz appliquée par 
Heisenberg. Toutefois comme ce probléme d’intégration numé- 
rique est relativement simple puisqu’au lieu d’avoir affaire A un 
systeme on a affaire 4 une seule équation a une seule fonction in- 
connue, nous exposerons d’abord une méthode générale d’inté- 
gration des systémes d’équation intégro-différentielles. 


2. Méthode générale d’intégration des systémes d’équations 
intégro-différentielles (#8), — Fock (38) a développé une mé- 
thode de résolution numérique des équations du champ self- 
consistent par approximations successives, mais cet auteur con- 
sidére des équations simplifiées dans lesquelles les coefficients 
7izx Non diagonaux sont nuls. Le nombre des paramétres 4 déter- 
miner est ainsi réduit ce qui facilite Vintégration des équations, 
mais les fonctions obtenues n’étant pas en général orthogonales, 
on ne dispose plus de paramétres ajustables pour les orthogona- 
liser rigoureusement tout en satisfaisant aux équations. On est 
done obligé de les orthogonaliser en cherchant A minimiser leur 
désaccord avec les équations différentielles, et le processus d’ap- 
proximations s’en trouve évidemment ralenti. Hartree (#4) a 
remarqué ce défaut de la méthode de Fock et n’a pas supprimé 
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les 24, non diagonaux, par contre les procédés de tatonnements 
et d’essais quwilemploie méritent & peine le nom de méthode. En- 
fin ces deux auteurs paraissent ignorer la méthode de Milne (°) 
pour trouver les fonctions propres et les valeurs propres d’une 
équation de Schrédinger. Nous allons donc exposer la méthode de 
Fock tout en conservant les 2, non diagonaux et en utilisant la 
méthode de Milne quand ce sera nécessaire. 

Nous supposerons tout d’abord ’hypothése (62) réalisée, ce qui 
réduit A une seule : r, les variables indépendantes. 

Le systéme a intégrer est le suivant : 


ad? i ! —/ 
(65) Ci at + Qu(r)fi= fi + ~Y hiefk —»>' Qix(r) fs 
k k 


dans lequel l’accentuation du signe = signifie que la sommation 
doit étre effectuée pour toutes les valeurs de Pindice k sauf k = t. 
Les fonctions Qy (r) ont pour valeur 


(66) Qn = [ fi*(r') Gaalr, r')fe(r’) - dr’. 


Les noyaux Gj (r, r’) dépendent des lois de force considérées. 
La fonction Qi (r) se compose d’un terme analogue a Qiz (7), et 


d’un terme 7 qui provient du développement du laplacien 


en coordonnées sphériques. 
Les fonctions fi solutions du systéme (65) doivent étre orthogo- 
nales et normalisées. 


+a 
0 fi*(r)- felr)-dr = Fie. 


A. Formation d’une premiere approximation. — Pour former 
une premiére approximation on utilise la méthode de Ritz comme 
Va fait Heisenberg : 

a) On prend comme fonction d’onde partielles de premiére ap- 
proximation les fonctions propres d’un oscillateur harmonique de 
fréquence arbitraire. 

b) Avec ces fonctions d’onde on forme les densités mixtes et 
Vénergie E du noyau. 

c) On détermine la fréquence « de Voscillateur de fagon a mi- 
nimiser |’énergie. 


oe 
ie) 


YE 
aa” 


| 


== 0. 


Y 
R 
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B. Utilisation de cette premiére approximation. — Nous désigne- 
rons par fi(r) les fonctions approximatives connues. En POrAnG 
ces fonctions dans (66) on calcule les Qi (r) qui deviennent ainsi 
des fonctions connues que nous désignerons par Qy(r). | 

Le systéme (65) s’écrit alors : 


(67) CFE + Welt — i= Dah —D Taf 
k k 


Enfin nous porterons également les fonctions de premiére ap- 
proximation dans le membre de droite de (67) de fagon a trans- 
former ce systéme d’équations différentielles simultanées en un 
systéme d’équations indépendantes. 


ad? Sle ie sae = — 
(68) TE + Qulr)fi — afi = Lah —Y Gat 
C. Intégration du systéme (68). — Nous allons intégrer séparé- 
ment les équations qui forment le systeme (68). Pour cela posons 
dg = 49 + 231 


et faisons passer le terme ),! f: dans le membre de droite en rem- 
placgant la fonction inconnue fi par la fonction de premiére appro- 
ximation fi. 

On obtient ainsi 


Tre 
(69) C1 + Wile) fi — 2th = wa) 
avec 


wir) = fi + Yueh —Y Gain 
k k 


L’équation (69) n’a de solutions satisfaisant aux conditions li- 
mites que si le membre de droite est orthogonal aux fonctions 
propres du membre de gauche égalé a zéro. Sid est valeur propre 
de l’équation homogéne obtenue en égalant 4 zéro le membre de 
gauche de (69) 


hi ta 
(70) O58 + Cal — nf = 0 


cette équation n’a qu’une fonction propre que nous désignerons 
par f° et il faut que : 


+o 
(71) [ ies deat 


a 
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On pourra toujours satisfaire 4 cette relation en ajustant conve- 
nablement la constante ii! contenue dans »;(7). 

Pour trouver la valeur. propre :° de l’équation homogéne (70) 
et la fonction propre correspondante fi? le procédé le plus str 
consiste & appliquer la méthode de Milne qui est exposée plus 
loin. La fonction f° peut étre normalisée sans ditficulté puisqu elle 
n’est définie qu’a un facteur constant pres. 

Connaissant la fonction fi on peut, suivant la méthode de Fock, 
intégrer l’équation non homogéne (69) en utilisant les fonctions 
de Green généralisées. 

Pour abréger, posons 


Qii(r) —_— iP 1 bo 
L’équation (69) prend la forme 
afi 


(72) Ci xi Pfi = i(r). 
Considérons |’équation 

dF 
(73) Ci Fa a PF = f?. 


Soit F; (r) une solution de cette équation (73) nulle pour r = 0 et 
quelconque pour r = ~. 

Soit Hi(r) une solution de (73) nulle A Vinfini et quelconque a 
Vorigine. On prendra toutefois ces fonctions Fi (r) et Hi(r) dis- 
tinctes. | 

Remarquons que ces fonctions sont faciles & déterminer par In- 
tégration numérique. Au moyen de ces fonctions Fi(r) et Hi(r) 
on construit une fonction de Green de la fagon suivante : 

\ G(r, 7’) = Fir) f(r’) + Hil’) fP(r) pour r< re 
) G(r, r’) = Far’) - f(r) + Hil) for’) pour r<r 
Cette fonction est continue en tout point sauf en: r = ral 

En utilisant cette fonction de Green la solution particuliére de 

’équation non homogéne (72) se met sous la forme 


+o 


(75) gir) =— {8 G(r, r’)-wa(r’) dr’ 


(74) 


r +a 
= — Hi(r) f foo(r’)-wa(r’)dr’ — Filr)> J fio(r’) -wi(r’)ar’ 


( Jo 
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Pour établir que gi(r) est bien solution de l’équation (72) 
remarquons tout d’abord que 


d[.dy  dfe dy | 
(76) ; C5, |G —y S| = 14S + Pylf® 


est une identité quelle que soit la fonction y. 

Pour y = Fi(r), ou y = Hi(r) le second membre a pour va- 
leur : [f°}*, et par conséquent en intégrant les deux membres de 
(76) et en tenant compte des valeurs aux limites de Fi(r), Gi(r) 
et fi®(r) on obtient : 


| C, [ios hae a = [trea 


dH, df,® +e 
lo [re ini — HF] =— F [fio }’dr 


(77) 


I] suffit maintenant de différentier deux fois Pexpression (75) 
et de simplifier le membre de droite en tenant compte de (71), (77) 
et de la normalité de fi° (r). 

On trouve 

dg; 
Ose + Pgi = wi, 


done gi(r) est bien une solution particuliére de (72) et la solution 
générale a la forme 


fir) = gir) + afi(r). 


D. Orthogonalisation et normalisation des solutions du systéme 
(68). — w(r) étant une fonction linéaire. des paramétres ii) et 
Ak, et d’autre part la fonction gi(r) se déduisant de w:(r) par une 
opération linéaire, on voit que les fonctions fi(r) sont linéaires en 
a, di, du, et ont la forme : 


fir) = af %(r) + afar) + SD uafit(r). 
k 


Chacune de ces fonctions doit satisfaire aux conditions sui- 
vantes : 


n+ 


J, 


ee 
; fi*(r)- filr)dr = bin, 


” fiir) Bar = 4 
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et d’autre part les paramétres : «, 2s, 24, sont liés pour chaque i 
par la condition 


? 


*+4 00 
if wir) fi(r)dr = 0. 


On voit facilement qu’il y a exactement autant de conditions a 
remplir que de paramétres A déterminer, et par conséquent les 
fonctions fi(r) sont rigoureusement définies. 


E. Suite du processus d’approximations successives. — Pour 
continuer l’intégration du systéme il faut porter les fonctions fi(r) 
qu’on vient de trouver dans les intégrales Qi(r) et Qi(r), et dans 
tout wi(r). Ces opérations représentent d’ailleurs la partie la 
plus fastidieuse du calcul numérique. 

I] faut ensuite intégrer le systeéme (69) ainsi formé, et recommen- 
cer la suite de ces opérations jusqu’a ce que la différence entre les 
fonctions d’onde initiales et les fonctions d’onde finales soit infé- 
rieure a l’erreur qu’on se permet. 

Remarquons avec Fock que bien que le coefficient Qii(r) 
de ’équation homogéne (70) chanye a chaque approximation il 
n’est pas nécessaire de résoudre cette équation et de déterminer 
Ai, fi(r), Fi (r), et Hi(r) de nouveau chaque fois. Il suffit de déter- 
miner ces paramétres et fonctions une fois pour P approximation 
initiale, et ensuite d’inclure dans le cété droit de Péquation in- 
homogéne la différence entre la valeur initiale et les valeurs sui- 
vantes de Qi(r). 


3. Transformation du systéme d’équations intégro-différen- 
tielles en un systéme d’équations intégrales (°°), (37). — I] est inté- 
ressant de remarquer que le systeme d’équations intégro-différen- 
tielles (65) peut étre transformé en un systéme d’équations inté- 
grales par l'utilisation d’un procédé analogue a celui qui nous a 
servi pour intégrer (69). 

Pour cela nous séparerons dans Qi(7) la partie intégrale, que 

li(li + 1) 
2 


nous désignerons par le méme symbole Qu(r), du terme ieags FPS 


provenant du développement du laplacicn. 
En divisant par Ci et en posant 


Mi 
Coie 


M. Marricon 5 
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Je systéme (65) s’écrit : 


@ (ls + 4) | 
(78) oe — ei — “Ah = wir, 
avec 
h Qi 
gi(r) = c.f as fe. 
k k 
Posons : 


e = kr, 
et considérons l’équation homogéne obtenue en égalant a zéro le 
membre de gauche de (78) 


d*fi li(li + 1) 
(79) ea hk— EP pao. 
Cette équation a une solution Fi(e) nulle & V’origine : 
eg ANAS 
Fi,(e) = git( =) (o-1she) 


ott 02 

"WO = Ta te HY + aan Fa) th 
et une solution nulle & l’infini 

Gye) = (— Ayyghtt( 2 Mg teP) a oF, 
Ces deux fonctions satisfont a: 
dF. 
a ou 

Avec Fn(p) et Gu(e) on forme une fonction de Green qui per- 

met d’écrire la solution unique de Péquation (20) sous la forme 


(80) 


(81) fi(r) — —{ F,,(kir’) . G,(kir) -gicr’)dr’ 
iis | F;, (ker) Gy,(her’)«gi(r’) «dr’. 


On vérifie facilement que (81) est bien la solution de (78) en diffé- 
rentiant deux fois et en tenant compte de (80) et du fait que Fy 
et Gi satisfont a l’équation homogéne (79). i 

Ainsi le systéme d’équations intégro-différentielles (65) est 
remplacé par un systéme d’équations intégrales non linéaires (84). 


4. Intégration de ’équation du noyau d@hélium. — A. Unités 
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et interactions employées (2%), (**). — Pour effectuer les calculs 
numériques nous emploierons les unités de Feenberg (?%). 


; h? 
unité de longueur \/ intima, = 8,97-10- cm. 
unité d’énergie : mc? = 8,119-10-7 ergs. 
D’aprés Inglis les lois de force sont : 


Gat = — 5 Co 1 Fe 
eet EES Peay — b? | r— ra 
G,*? = 0 





Gp = — Por mir 
G,? = 0 

G,t = 

pi=0 

Jt = Bil — ge P 
js = 0 


J4 = Bge BP 
les coefficients ont les valeurs suivantes 
Bea 72 === 435 
g=0,205 b= 16, 
et le terme de Coulomb a pour valeur : 
0,312. 
"fe 
B. Premiére approximation. — Pour obtenir une premiére ap- 
proximation nous utiliserons la méthode de variation de Ritr. 
Nous trouverons en méme temps la valeur de premiére approxi- 
mation de l’énergie 4 laquelle nous voulons comparer la valeur 
finale. 
Dans le noyau d’Hélium 
Pi = Po =m = Ng = 1. 
I] n’y a done qu’une seule sorte de densité 
(|p| x) = $*(z)-¥-(z’). 
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Nous négligerons la différence entre la masse du proton et celle 
du neutron, ce qui donne a !’énergie (57) la valeur 


Li penis 
ah Y* (a!) = Y(ar')» = op *(x) + Yr) dander 
sea! r | 


b2 aed ny Zi 2 . 
+ 2G + BE— g)} { ¥*()-v(a) se PIP? Pys(ay -a(e dade’ 
Prenons comme pcint de départ, la fonction d’onde d’ordre 
zéro de l’oscillateur harmonique : 


Wr) = be 2”, 


La condition de normalisation 


24 co 
a) V(r) -4nr2dr = 1, 


fixe la valeur du coefficient 4, et on obtient : 


(82) Wr) = (S\eaae 


avec cette fonction d’onde nous allons calculer les différents 
termes de E. 
1° Energie cinétique 


We (PB 2 
Ecin = 3H, = 3 [ w) <— gan (Fa +5. i) Yr) Sertdr, 


Avec nos unités : 





h2 Ot 4 
Sx2m”— 2 
On obtient 
Kein = 2,254, 


2° Energie électrostatique des protons 


Hea [ vr)%dr- sin 9-d0-dp. [ 0) -r’?dr' sin 6'd0' «dy. 


V r? r'2__Ipr' cosh’ 


Ey = 22 / 5 = 0,252 
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3° Energie d’échange entre protons et neutrons 


wr 
Eon = AC. + B(2 — lj Y(r)-r? sin 6-drd§-dy- [ve 
e— AU + r'2—2rr! cos B)p2gp'. sin 0'd0'dy’. 
at ac + B(2 — g)]2 2 
(a + 2b2)'!2 
Nous obtenons donc la valeur suivante de l’énergie de liaison 
343,482'!2 
(2 + 32)°/2 
On minimise cette expression en annulant sa dérivée par rapport 
az 


BE = 2.252 + 0,252 — 


dE ) Bl 
a, = 2425 + eG 8] ie O 
V2 (x + 32)/2 





La racine de cette équation qui correspond 4 un minimum pour 
E est : 


a = 416 
La valeur de l’énergie qui lui correspond est : 
E = — 50,75 


Elle se décompose de la fagon suivante : 
| Energie cinétique : 93,6 
(83) Energie proton-proton: 1,61 
Energie proton-neutron : — 145,96 


C. Deuxiéme approximation. — Nous allons intégrer numéri- 
quement maintenant par approximations successives, l’équation 
du noyau d’hélium. 


Nous traiterons l’action de Coulomb, comme une perturbation 
et nous négligerons son influence dans l’équation d’onde, ce qui est 
légitime d’aprés les résultats (83). 

Nous négligerons également la différence entre la masse du 
neutron et celle du proton. Dans ces conditions 

H! = H?= H 
et ’équation du noyau d’hélium s’écrit : 


(84) 2HY(x)+-2[C+B(2—8)] i ye(a') Vat) ent Pada! -Ya)=M4(@) 
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4 
ou en remplagant H par sa valeur : —% + 3 A et en supposant 


que | ne dépend que de:r 
Ay + H(r)-4 = py 


avec : 


et: 


H(r)= Mike 3 X 2[C-+B( (2— oi f (r’)- EP (A+ r2—2rr! cos 8!) padp! . sin6’d9'd¢’. 


+ © , ' 
“aif V(r’ fe lg PEA dr", 
avec : 


a= —2 x JC-+B2— g)]. 


Comme d’autre part : 


a? 2d 
Dan ey 


L’équation d’onde prend la forme définitive 


(85) (aa + 2 %)e0) + HO)-90) — wr) = 


Nous devons chercher une solution Y(r) de l’équation (85) qui 
soit partout continue, finie, univoque et nulle a l’infini. Cette so- 
lution n’existe que pour certaines valeurs de » appelées valeurs 
propres que nous chercherons également. 


Pour faciliter ’intégration de pepe hie (85) faisons le change- 
ment de fonction : 


o(r) = r-Y(r) 
Ona: 


dy & “ta 


de? —"lae +7 a 
Multiplions (85) par : r, nous obtenons : 
a2 
(86) a + Hire — np =0. 


La solution de cette équation doit étre nulle pourr = + o, mais 
aussi pour: r = 0. 
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L’équation (86) est seulement d’apparence linéaire, car la fone- 
tion : H(r) 


+o 
(97) = LM yer pte tee — ete a 
0 


contient la fonction inconnue (7). 

Pour intégrer l’équation intégro-différentielle (86) nous suivrons 
le processus d’approximations successives habituellement utilisé 
dans la méthode du champ self-consistent : nous formerons la - 
premiére fonction H,(7) en introduisant dans l’intégrale (87) la 
fonction d’onde de premiére approximation (82). L’équation (86) 
ne sera plus alors qu’une simple équation de Schrédinger dont la 
résolution nous donnera une fonction d’onde de deuxiéme appro- 
ximation. Cette fonction d’onde sera portée a son tour dans (87) 
et Ja fonction H,(r) ainsi obtenue permettra la formation d’une 
équation de Schrodinger dont Vintégration fournira la fonction 
donde de troisieme approximation, et ainsi de suite jusqu’au 
moment ou la fonction d’onde de pee approximation coinci- 
dera avec le nombre de chiffres significatifs voulus avec la fone- 
tion d’onde de l’approximation précédente. 

A chaque stade l’équation (86) est une équation de Schrédin- 
ger qui, étant homogéne, définit la fonction g a une constante 
multiplicative prés. Cette constante sera choisie de fagon & satis- 
faire a la condition de normalité. 


+a +o 
| | vir) pdartdr = 4 [ | ¢(r) |2dr = 1. 
0 0 


En portant la fonction (82) dans Vintégrale (87) on obtient 


ab? 
ax’ |s ee A HE 


(a+ byl” 
= 281 091156", 


‘a eee 


Nous devons donc intégrer Véquation de Schrédinger 
Py J 09 eS OE tee 1) mae eh 
(88) ar? + 281,09 e spe LT 6 : 
Pour cela nous emploierons la méthode de Milne. 


D. Méthode de Milne (°) pour lintégration numérique et la déter- 
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mination numérique des valeurs propres d’une équation de Schro- 
dinger. — Soit l’équation 


d* 
(89) aT t Ge, *) =0, 


qui contient comme cas particulier l’équation d’onde du type 
8x2 


du 
1 + 7 [W — V(x)]u = 0. 


Supposons que l’intervalle soit —«o <a4< + o. Soit a, une 
valeur de x dans ,l’intervalle, et soient u,(x) et u,(”) deux solutions 
de l’équation (89) satisfaisant aux conditions 


U(X) = 1 Ug(X) = 0 
Uy'(%9) = 0 Us (Xo) = a 0. 


La constante a@ est arbitraire et dans chaque cas particulier 
doit étre déterminée de facon a rendre le travail numérique aussi 
simple que possible. Ces solutions satisfont a Pidentité bien connue 


Uy(X) + Ug’ (x) — Ug(x)- Uy'(x) = a. 
car d’une part 


Uy (Xo) *Uy'(Xp) — Ug(Xp) - Uy'(%) = a 


et d’autre part : 


Pe [Wy - Wy! — Ug: Uy'] = Uy- Uy" — uy- ty". 
= Uy,-U,(G(x, 4) — G(x, d)) = 0 
Soit une fonction (x) définie par 
w(x) = [2,2(x) + u,%(a) 12 
Différentions cette équation deux fois 


w'(2) = Uy°Uy + Uy* Uy’ 


Ww 
w"(2) = wy? + w + w'g? + Wyuy" + Ugus! (UyUy' + UgUy')® 


Ww ws 
SSH) G(x, )(u,? + u,%) 4 (a's? + u',?) — (UyUy' + UgU,')? 
a NE a a) ty = Ugtts )* 
W ws 


2 
= — G(x, yw +S. 


Donec (zx) satisfait a l’équation 


d? 2 
a + Gla, iw —=0 
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Les conditions initiales sont : 
W(X») = 1 
W'(%) = 0. 


Maintenant la solution générale de (89) peut étre exprimée sous 
la forme 


%% 


u(x) = Cw(2) - sin } a [wae ——i te ! 


dans laquelle c et « sont les deux constantes arbitraires nécessaires, 
dont l’une C est multiplicative comme il se doit puisque l’équation 
(89) est homogéne. 

En effet, on a 


du dw, 
Betas ants UES ein fF GCE « COBis. f+ | 
d*u d*w w 
— — (| ——.. gin |... ! ae AS erah SACD 
ie Cr sin } Pag Dee cos | (aw 
dw : 
— aCw-? — cos |--- ‘ — Caw sin |---| 
dx 
ee ee 
== (sini 3st xX “ro |- 
et par conséquent 
du : ; d2w a 
de + Gl iu = Gsin }---(x|F3 + Gle, Nw —% , 


==. 0 d’aprés (106). 


Puisque avec les hypothéses faites w(x) ne s’annule 4 aucune 
des deux extrémités de lintervalle, il est clair qu’une solution 
satisfaisant aux conditions désirées aux deux extrémités de l’in- 
tervalle ne peut-étre trouvée que si 


+a 
(90) - th w4dre =n 


T ]—o 


dans laquelle n est un entier positif, car 
+2 
ut+o)=0 si { w dz —a = Nyt 
X% 
u(i—_«o) = 0 si i dz — a = nr 
a) 


et en éliminant simultanément 2, et « entre ces deux relations on 
obtient (90). 
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Maintenant |’intégrale 


+0 
N = a w dx 
beh ipa 


est une fonction croissante de i, telle que si nous choisissons k 
valeurs de 2, soient A, < ?, << --- <2, nous obtenons k valeurs de 
N: N,<N,< ---< Ny. Si un entier n tombe dans l’intervalle 
N, << Ny, nous pouvons utiliser la méthode d’interpolation 
pour obtenir approximativement la valeur correspondante de ), 
A = dn qui sera précisément la ni*™ valeur caractéristique de ) 
comptée par ordre de grandeur. 


E. Intégration de l’équation (88). 


2, : 
(91) Te + 281,09¢- 11.55%» ug 9, 


L'intervalle étant 0 et + 0, nous appliquerons la méthode de 
Milne en cherchant a rendre l’intégrale 


Pr 
- f w-*dx égale a l’unité. 
/0 


Nous prendrons : ry = 0, a = \/n. 
Nous devons done intégrer numériquement Péquation 


dw ak Bp 
apt = (v — 281,09e— 11,56") q, 4 = 


avec les conditions initiales 


pour certaines valeurs de p. 
Pour obtenir une valeur de p. satisfaisante rappelons-nous que 


ms te 


3 


Pour trouver la signification de 2, multiplions a gauche (84) 
par Y* et intégrons, on obtient : 


NADAS | 
A= 5 Eein of Eon 
soit, en utilisant les résultats (83) 


h=— 99,16 
» = 132,21. 
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Comme notre calcul doit améliorer la valeur de l’énergie nous 
essaierons 
w= 192.5. 
Intégration de lV équation 





132,5 
3 


w" = 132,50 4+ — 


— H(r)w. 
Les conditions initiales sont 

w(o)=1 

w'(o) = 0 

Nous prendrons la suite des valeurs de r équidistantes de 
Ar = 0,005. 

La méthode la plus simple pour intégrer une équation différen- 
tielle est la méthode de Milne. 

Méthode de Milne ('), (7). 

Les quatre premiéres valeurs de w, w’, w’’ étant obtenues au 
moyen de quelque méthode appropriée, telle que celle de Runge- 
Kutta, les valeurs sivantes sont trovvées ainsi : 

1° On trouve la premiére approximation de la valeur siuvante 
de w’ au moyen de la formule 


4Ar "1 uf ” 
Wn'D = W'n_4 + 7 (20 Nils ee Ne a 2w na): 
2° La valeur suivante de w s’obtient au moyen de 
Ar / / / 
Wn = Wn-g + = (Wn'® + 40'n_y + n-2)+ 


3° En portant cette valeur de w dans ]’équation a intégrer 


132,5 


wy 


w" = [132,5 — H(r)]w + 





on obtient w,’’. 
4° On trouve la deuxiéme approximation de w’, en substituant 
cette valeur de w,'’ dans 


’ Ar " / ” 
Wn = w n-g + in (Wn + 4w'n_y +H n-2)- 


Si wn’ et wa’® coincident au nombre désiré de chiffres signi- 
ficatifs, nous pouvons considérer w,'® comme correct, et recom- 
mencer le méme processus pour l’intervalle suivant. 

Si wa’ et wa’® ne se correspondent pas trés bien et qu’aucune 
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erreur ne puisse étre trouvée dans les calculs, alors calculer la quan- 
tité 
| @n'@ — wy'® | 
Th no a ae 
qui est l’erreur due a la méthode. Si cette quantité est assez grande 
pour affecter le dernier chiffre significatif que nous voulons ob- 
tenir, la seule chose a faire est de diminuer l’intervalle Ar. 

La méthode de Milne est trés simple, et 4 chaque étape l’examen 
de la valeur de < sert de contréle précieux pour l’exactitude des 
calculs. 

Méthode de Runge-Kutta (7). — Nous appliquerons la méthode 
de Runge-Kutta pour trouver les quatre premiéres valeurs de w 
et de w’ nécessaire au démarrage de la méthode de Milne. 


Dans la méthode de Runge-Kutta les accroissements Aw et Aw’ 
sont calcnlés au moyen d’une suite déterminée de formules. 





w = y' u' = (132.5 ae H(r))w + = — f(r, wv) 
k, = f(r, w)Ar l,=u-Ar 
he f(r +3, w+a)ar = (w+ )ar 
A I | k 
b= f(r+ Ze 4B)ar b= (u+2)ar 
ky= f(r + dr, w + 1,)Ar ly = (w+ kg) Ar 


1 
Au= (ky + Qhy+ ky th, Aw= a(t + Ql, + 21, + 1,). 


Comme nous avons calculé les valeurs de H(r) pour des valeurs 
de r différant de 0,005 la méthode de Runge-Kutta nous donnera 
les valeurs de w et w’ pour r = 0 — 0,01 — 0,02 — 0,03, et il est 
plus simple de calculer les valeurs intermédiaires avec une for- 
mule dinterpolation, telle que la formule de Bessel qui utilise les 
différences centrales. 


Ayo 1 A*y_, + Aty, 9 Agy_2 + Aty_, 
8 2 


Yos = Yo + 9 3 + 798 a ae 


L’intégration terminée, nous calculons 


132,5 (70:80 
EE Ra 


avec les colonnes 14 et 12. 
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Cette intégrale a pour valeur: 1,000.464. Selon la méthode 
de Milne nous devons recommencer toute Vintégration avec 
» = 133, ce qui donne une vaieur de I inférieure a 1, et par inter- 
polation on déduit la valeur de » qui est valeur propre. 

Connaissant, par interpolation également, la fonction w qui cor- 
respond a la valeur propre » = 132,81 on en déduit la fonction 
propre non normalisée, solution de (91) en calculant : 

r 
9, = w(r) + sin jaa w 2dr ‘ 


La fonction d’onde ¥,(r) correspondante s’obtient d’apres (97) en 
faisant : 


,(7) ee 


La valeur initiale 1(0) qui est indéterminée s’obtient par extra- 
polation des valeurs suivantes. 


1 
Enfin le calcul de if 4 |9,|? dr permet de normaliser (. 
0 


F. Troisiéme approximation. — Pour trouver une troisiéme 
approximation nous devons recommencer la méme suite de cal- 
culs que pour la deuxiéme approximation, mais en partant de la 
fonction d’onde ) que nous venons d’obtenir. 

La principale difficulte réside dans le calcul de 


ie ‘\2 ‘2 
H(r) = ele [ | v(r’) ign shy e— 16(7 + r’) yr -dr’. 


Cette fonction doit étre connue 4 10° pres pour toutes les va- 
leurs de r. Il faut donc calculer l'intégrale en prenant des valeurs 
de r’ suffisamment rapprochées, et cela pour chaque valeur de r. 

Pour raccourcir ces calculs nous désignerons par ,(r) la fone- 
tion donde de premiére approximation, par ,(r) la fonction 
d’onde de deuxiéme approximation que nous venons d’obtenir. 
Posons 

a[y7] = [42 — (hI? 
H,(r) se décompose de la fagon suivante 
H,(r) = H,(r) + 6H,(r), 
avec 
Hy) = one apee(r'y) fe tor Ph 88 +1) 


Comme 5H,(r) est del’ordre du 1 /100 de H,(r),on voit qu’on obtien- 
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dra H,(r) avec l’exactitude voulue en calculant oH rye Os 

Pour obtenir ce résultat, il suffit de prendre des valeurs de r’ 
plus espacées, et des valeurs de r également plus espacées, les va- 
leurs intermédiaires de 3H,(r) se déterminant facilement par in- 


terpolation. 
Pratiquement pour calculer 3H,(r) on décompose ainsi : 


Hy) = 2 


G(r) = I,(r) — I(r) 
+ 
I,(r) = 89,92 if r’-[d(r') Pe 16")? gpr 


+00 
I,(r) = 20,02 [ 1’ -8[U(r')]2-e— 16(" +) ap 
0 


Connaissant H,(r) = H,(r) + 3H,(r), on a les éléments né- 
cessaires pour calculer la troisiéme approximation. 
Apres quelques essais, qui donnent des valeurs de I trop faibles, 
on intégre l’équation 
Pw : 
da = [&— Hr) + 4, 
pour 
w= 4136 
et 
w= 137, 


ce qui fournit deux fonctions w, et ,. 


On calcule 
136 (+e 
I, ae Vv o if W,2dr 
0 





Tv 


137 Bi 
ay 


/ 





I, —s w,"dr, 


En adoptant la formule d’interpolation 


Da Viv(a— b,), 


on trouve que I = 1 pour: BU = 136,25. 

Une simple interpolation linéaire permet de calculer W136,95 CON- 
Nalssant W436 et W5.. 

A partir de cette fonction ®, on trouve / comme pour la pré- 
cédente approximation. 
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Si on compare cette fonction d’onde a celle de la précédente 
approximation on trouve qu’au maximum 


1 MUR IN 
y 1000 
or cette variation relative de | est due a dH,(r). 

Si nous calculons 3H,(r) en suivant la méthode déja indiquée 
on constate que 0H,(r) est une fonction der qui, 4 un facteur cons- 
tant prés, est identique 4 3H,(r). 

Ce facteur constant est tel que : 





Par conséquent, une nouvelle intégration modifierait ) au maxi- 
mum de 4.10~*, ce qui est une précision supérieure A celle qui est 
nécessaire pour apprécier la méthode du champ self-consistent. 

On trouve de méme que ». serait modifié de + 0,29. Nous 
adopterons done la valeur : »p = 136,54. 

Nous pouvons donc nous arréter 4 la troisiéme approximation. 

Dans le tableauI nous donnons les valeurs de la fonction d’onde 
initiale obtenue par la méthode de Ritz, ainsi que notre fonction 
d’onde finale. La figure 1 représente les courbes représentatives 
de ces fonctions. 


G. Calcul de ’énergie. — Pour déterminer ]’énergie nous avons 
dune part : 
3 1 
A= sen == 5 Bein + Epn. 

Donc : 

; Kein + Epn = : 136,54 

,) cin + 8 Soe 4 x ’ 

= — 102,405 


et d’autre part : 


Ey = — | ¥*(r)+H(r)-4(r) -4er2dr 


= — f drt-43(%)- Hip) 
= — 149,90. 


Par ailleurs nous ne ferons pas d’erreur sensible en admettant que 
l’énergie proton-proton qui est en elle-méme tres faible, a varié 


TABLEAU I 


Y initial y initial 


6 , 882 
6,814 
6, 
6,588 
6,439 
6,267 
6,074 
5 , 864 
5,637 
5,3 


mooooeooooo S&S 
COONAUP ON = 


~ «=~ we we «0 (ee, & 


SDD DH HD DH OVO OU OF O71 ie OO bt 
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0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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dans la méme proportion que l’énergie proton-neutron. Au total 
nous obtenons done pour l’énergie du noyau d’hélium 


= — 53,26 
au lieu de : — 50,75 par la méthode de Ritz. 
L’écart par rapport a la valeur exacte : — 55 a done diminué 
gern. 
L’énergie du noyau se décompose de la fagon suivante 
Eun = 94,99 
Ee = 4,65 
Epn = — 149,90. 


N. B. Dans un précédent calcul nous avons intégré Véquation 
du noyau d’Hélium dans le cas ov il n’existe qu’une interaction 
du type de Majorana entre neutrons et protons, et Vaction de 
Coulomb entre protons. Les coefficients de la loi de force étaient 
ceux donnés par Feenberg (”). 


a=174 
b? = 48. 


La valeur de l’énergie s’éloignait beaucoup plus de la valeur 
réelle que dans le cas actuel. 
On trouvait : " 
E = — 34,97 
par l’intégration numérique au Jieu de : 
E = — 31,25 


par la méthode de Ritz. 
La modification de la fonction d’onde était presque identique 
a celle que nous venons de trouver et qui est représentee figure 4. 


5. Valeur de la méthode du champ self-consistent dans l'étude 
des noyaux atomiques. — La différence entre la valeur calculée 
de l’énergie du noyau d’hélium et la valeur réelle — 55 montre 
que la méthode du champ self-consistent n’est pas rigoureusement 
applicable aux noyaux atomiques. 

Toutefois, il convient de remarquer que énergie calculée est 
la différence de deux termes assez grands : 94,99 eb — 149,40. 
Par conséquent un léger changement dans la loi de force reta- 
blirait ’énergie A sa valeur réelie. Pour regagner les unités mc? qui 


M. MATRICON 6 


82 APPLICATION DE LA METHODE DU CHAMP SELF-CONSISTENT 


manquent il suffit d’augmenter l’énergie potentielle de 1,16 .%. On 
peut y arriver trés facilement en augmentant par exemple les 
coefficients B et C de la loi de force. 

Il est done raisonnable d’admettre que la méthode du champ 
self-consistent doit fournir une description assez juste de la réalité 
physique puisqu’il suffit d’une augmentation de 1,16 % de la valeur 
d’un paramétre multiplicatif pour que les résultats de cette mé- 
thode concordent avec les résultats expérimentaux. 

Il parait par conséquent possible de demander a cette méthode 
des renseignements qualitatifs sur la structure des noyaux légers. 

Il est toutefois nécessaire de rechercher un mode de calcul sus- 
ceptible de donner des résultats quantitativement bons. 

La méthode qui se présente immédiatement a l’esprit est la mé- 
thode de perturbation telle qu’elle est exposée page 42. 

Pour que les calculs se fassent simplement ii y a lieu de Pappli- 
quer au probléme du champ self-consistent résolu par la méthode 
de Ritz. C’est ce qu’a fait Inglis (22) au sujet des noyaux d’Hé- 
lium et de Lithium. 

Dans la méthode de Ritz-Heisenberg les fonctions d’onde indi- 
viduelles des particules nucléaires sont les fonctions propres de 
Voscillateur harmonique dont on détermine la fréquence pour 
minimiser l’énergie. 

L’Hamiltonien qui correspond a une telle solution est done de 
la forme 


N 


> [5 Ricks ire ne 


i=1 


Le paramétre k étant relié trés simplement a la fréquence de l’os- 
cillateur harmonique. 

C’est la différence entre cet Hamiltonien et ? Hamiltonien réel 
du probléme nucléaire que l’on doit introduire comme perturba- 
tion du probléme self-consistent. Le restant du calcul se conduit 
suivant la méthode générale indiquée pages 42 et suivantes. 

I] semblerait, d’aprés les derniers resultats, que la méthode de 
perturbation appliquée de cette fagon, ne converge pas rapide- 
ment dans le cas des noyaux plus lourds. I] parait donc intéressant 
d’appliquer dans ce cas la méthode du champ self-con sistent 
sous sa forme exacte, ce qui améliorerait déja les résultats, et de 
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conduire le calcul de perturbation 4 partir des fonctions d’onde 
ainsi trouvées. 
On pourrait simplifier considérablement le travail numérique, 


tout en conservant les : de la précision de la méthode en se 


bornant a la deuxiéme approximation. Dans ce cas les opérateurs 
intégraux ne seraient plus des fonctionnelles, mais des fonctions 
déterminées, une fois pour toutes a partir de la premiére approxi- 
mation, et les équations intégro-différentielles deviendraient de 
simples équations de Schrédinger. 

Ce travail a été exécuté sous la direction de M. le Professeur 
Francis Perrin. Je ne saurais assez le remercier et lui témoigner 
ma reconnaissance pour l’intérét qu’il a porte a ces recherches et 
la bienveillance avec laquelle il m’a toujours conseillé. 

J’adresse aussi mes vifs remerciements 4 la Caisse Nationale 
de la Recherche Scientifique pour l’aide importante qu’elle m’a 
accordée ainsi qu’a la Compagnie Générale de Radiologie qui m’a 
laissé toute la liberté nécessaire pendant les travaux préliminaires 
et. dont le Conseil d’Administration m’a fortement encouragé et 
aidé. 
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